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Задача о бесконечной пластине, нагруженной нормальной силой, движущейся 


по сложной траектории” 


А. В. Галабурдин"” 


' Донской государственный технический университет, г. Ростов-на-Дону, Российская Федерация 
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А. У. СааБигат 


'Роп а Тесьиса1 Ошуегзйу, Возюу-оп-Ооп, Вазз1ап Еедеганоп 


Введение. Предлагается метод решения задачи о бесконеч- 
ной пластине, лежащей на упругом основании. На пластину 
действует периодическая нагрузка в виде силы, перемеща- 
ющейся по произвольной замкнутой траектории. Цель ис- 
следования — разработка численного метода решения задач 
теории упругости для тел, находящихся под действием по- 
движной нагрузки. 

Материалы и методы. Учитывая периодичность рассмат- 
риваемой нагрузки, она раскладывается в ряд Фурье на вре- 
менном отрезке, длина которого равна периоду нагрузки. 
Решение исходной задачи строится посредством суперпози- 
ции решений задач, соответствующих нагрузке, задаваемой 
слагаемыми описанного выше ряда Фурье. Окончательное 
решение задачи представляется в виде отрезка ряда. Каждое 
слагаемое при этом соответствует решению задачи о воз- 
действии на бесконечную пластину нагрузки, распределен- 
ной по замкнутой кривой (траектории движения силы). Для 
нахождения этих решений используется фундаментальное 
решение уравнения колебания бесконечной пластины, ле- 
жащей на упругом основании. 

Результаты исследования. Предложен новый метод реше- 
ния задач теории упругости для тел с нагрузкой, движущей- 
ся по замкнутой траектории произвольной формы. Решена 
задача о бесконечной плоскости, по которой с постоянной 
скоростью движется сосредоточенная сила. Определено, что 
траектория движения представляет собой гладкую замкну- 
тую кривую, состоящую из дуг окружностей. Рассмотрен 
характер изменения перемещений и напряжений вблизи 
движущейся силы. Изучено распространение энергии упру- 
гих волн. С этой целью выполнено вычисление координат 
вектора Умова — Пойтинга. Исследовано влияние скорости 
движения силы на длину вектора Умова — Пойтинга. 
Обсуждение и заключения. Метод применим и при рассмот- 
рении более сложных объектов (плиты сложной формы, 
слоистые плиты, вязкоупругие плиты). Его преимущество 
— экономичность, так как для построения решения исполь- 
зуются уже известные решения задач. Окончательное реше- 
ние выражается в удобном виде — как сумма криволиней- 
ных интегралов. Полученные результаты могут быть ис- 
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Маетаё апа Мешоа$. Слуеп Фе репо Ч1сцу оЁ Фе 1Тоа4 ипдег 
сопз14егайоп, и 15 десотрозе4 ш а Еоипег зепез ш а Ише 
и\егуа| \Позе [епо 15$ еаца| ю Ще 1оа4 рето4. ТБе зоайоп 
{о Ше опоша|! рго ет 1$ сопзкицщеа Бу зарегроз1оп оЁ Фе 
5о[аоп$ ю Фе ргоетз сотгезропашо {0 Фе 1оа4 зресНеа 
Бу Фе 1егт$ оЁ Фе Еоцпег зепез Чезсге4 абоуе. ТБе Йпа1| 
5оНоп о Фе ргоет 15 ргезеще4 а$ а зеотеп оРа зетез. ш 
{615 сазе, еасВ 1егт сотгезроп4$ ю Ше зо[айоп оЁ Ме ргоМет 
ог Фе ппрасе оп ап шНипе рае оЁа 1оа4 Ч1зитоще4 аопх а 
с1озе4 сигуе (фе {галестогу ог Ше Ююгсе тойолп). То Йпа Шезе 
5о[ийоп$, бе шп4датеща| зоаноп ю Фе едианоп о? уганоп 
оГап шйпие рае [ух оп ап @азис Базе 1$ изед. 

Кезеатгсй Кезий5. А пе\у тео 1$ ргорозеа г зо]уше ргоБ- 
1ет5 оп фе @азИсйу Феогу Юг Бо1ез у а 1оа4 ЮПо\злис а 
с1озе ра о{ агойгагу зВаре. ТБе рго ет оГап шйице р1апе 
а1оп> у\/в1сВ а сопсепгае4 Гогсе тоуез аё а сопз{апЕ зрее4 15 
5о[уе4. П 1$ а&етишеЯ ФфаЁ Фе Надесюгу оЁГ тойоп 15 а 
оо с1о5е4 сигуе соп$1$Нпе оЁ сисуаг агсз. Тре Бевау1ог 
ог 415р!асетеп{$ ап згез5ез пеаг а тоуше гсе 1$ соп$19- 
егеД. ТБе епегоу ргораганоп оГ ®е е|азйс \’ауез 15 заФеа. 
Еог 1$ ригрозе, е соог4таез$ ог Фе Отоу — Роупйпе уес- 
{ог аге са1сшаеа. ТБе еНес оЁ Ве Югсе тоНоп зрее4 оп е 
|епо оРфе Отоу — Роупйпх уесюг 15 шуезисаеч. 
Обсиуяопт апа Сопсияоп5. ТВе тео 1$ аррИса Ме \Пеп 
сопз14егие® тоге сотрйех обБ]есё5 (р1айез оЁ сошр]ех зВаре, 
[ауегеа р]а{ез, у1зсое|азис р!аез). [5 ауашаее 1$ рго а Шиу 
ушсе Фе Кпоуи ргоет зошйоп$ аге изе4 ю БиПА Фе зо№- 


Ноп. ТБе Нпа| Дес1з1оп 1$ ехргеззе4 ш а сопуешепе югт - аз 
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пользованы в процессе проектирования дорог. Изучение 
распространения энергии упругих волн от движущихся 
транспортных средств позволит оценить воздействие ука- 
занных волн на строения, расположенные вблизи дороги. С 
учетом данных о характере изменения перемещений и 
напряжений оценивается износ дорожного покрытия. 
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Введение. Исследование динамических явлений, обусловленных действием подвижной нагрузки, пред- 
ставляет собой актуальную задачу, имеющую прикладное значение (например, при решении вопросов развития 
транспорта). В работах, посвященных указанной проблеме, рассматривались различные задачи с подвижной 
нагрузкой. В частности, показано, как на полуплоскость или полупространство (упругое изотропное, трансвер- 
сальное изотропное, вязкоупругое) действует нагрузка, движущаяся по бесконечной прямой линии с постоян- 
ной скоростью. В этом случае при решении задачи вводится подвижная система координат, связанная с дви- 
жущейся силой, что позволяет исключить время из числа независимых переменных [1-5]. Часть работ посвя- 
щена изучению действия на бесконечную пластину или полосу (упругую или вязкоупругую) равномерно дви- 
жущейся по прямолинейной траектории нагрузки. При этом используется тот же метод исключения временной 
переменной или рассматривается квазистатическая постановка задачи [6-12]. Наибольший интерес представ- 
ляют задачи, в которых длина траектории движения нагрузки конечна, а сама траектория представляет собой 
кривую линию. Часто этом случае используется конечноэлементное моделирование движущейся нагрузки [1 |- 
13]. В ряде работ при решении подобных задач используются вариационные методы (в частности метод Рэлея 
— Ритца) [14—16] или разновидность метода Галеркина, позволяющая сводить задачу к обыкновенным диффе- 
ренциальным уравнениям. При этом рассматриваются разнообразные объекты приложения подвижной нагруз- 
ки (плиты, слоистые плиты, вязкоупругие плиты, полупространства — как изотропные, так и анизотроп- 
ные) [17—19]. В данной статье представлен метод, который развивает идеи, изложенные в работах [20—23]. 

Материалы и методы. Рассмотрим бесконечную пластину, лежащую на упругом винклеровском осно- 
вании, которая находится под действием нормально приложенной силы, перемещающейся по замкнутой траек- 
тории. 

Задача сводится к интегрированию уравнения движения пластины, лежащей на упругом винклеровском 
основании [14]: 


ты 
АИ +0 Е =—, (1) 

р 

ЕЙ? 
где И’ — прогиб пластины; ДО = 12—) ; Е — модуль Юнга; у — коэффициент Пуассона; й — толщина пла- 
—У 
> рй К 

стины; с “= г ; р — плотность материала; К = г ; К — коэффициент податливости упругого основания; Р — 


сосредоточенная сила, движущаяся по замкнутой кривой 7 с постоянной скоростью а. 
Введем координату 5, отсчитываемую от некоторой фиксированной точки кривой у. Тогда сила Р, дви- 
жущаяся вдоль кривой '/ со скоростью а, будет описываться соотношением Р=Р(5-—аё). Функция Р(5-—ай) 


рр 
периодическая по 1, с периодом Т = — ‚ где Ё — длина кривой у. 
а 


Решение. Рассмотрим установившийся режим. Разложим функцию Р(5-—а!) в ряд Фурье по перемен- 


ь [ 
ной {на отрезке | —-—; — |. В этом случае коэффициенты разложения предстанут виде: 


2а 2а 
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Г 
т 
акт, 
с: = | Ро -ае 1. 
_Г 
а 
Сделав замену переменной интегрирования в интеграле 5—аё = 2 , получим: 
Г, 
я 215 225 
7 е У. е У 
с. = | Р(еде а, 
[@ 


-7 


А, в 
где 4, — коэффициенты разложения в ряд Фурье на отрезке |= | : функции Р(2). 
Тогда подвижную нагрузку можно представить в виде ряда Фурье: 


1 — 2ип(5—а!) 
Ре-а) = У Че й | 
а 
К=—осо 


Учитывая линейность задачи, ее решение можно представить в виде: 


Я 


со 
И = А.И’, 5 (2) 
К=Ь—со 
где И’, — прогибы пластины, обусловленные действием вертикальной нагрузки, распределение которой вдоль 
К 


25 —Р ати 
кривой у описывается функцией е ^/, изменяющейся во времени по закону е Я 


Для определения И’, воспользуемся фундаментальным решением уравнения (1), которое соответствует 


кап 
т 





Р=5(х-ж)5(0?- уе" ‚ где о} = 


Используя принцип предельного поглощения и традиционные методы построения решений дифферен- 
д 
циальных уравнений, можно получить фундаментальное решение уравнения (1), которое при А > имеет 
б 
ВИД: 


и Ух, У) = [№ (%)- Ко (а>К) › 


1 2 | 
где К = |= } +(У- у у № = ^ 6, ; © = ве = о. ‚ Ко(2) — функция Макдональда. 
йе 


2 


© 
При А < —— решение имеет вид: 
с 


1 т 
= — НН (ув)- К, (78 


2 
— -К, Н ры (УК) — функция Ганкеля. 





- 21п(5—аг) 
Тогда И’, = фи (х, у, ху (5), У (5))е и. 45. Используя известные формулы теории тонких пластин и 


у 
полученную по приведенным выше соотношениям формулу, определяющую прогиб И’ (2), можно вычислить 


перемещения и, и, и напряжения с., о, и ох, в любой точке пластины. 
При больших А приходится вычислять интеграл от быстро осциллирующих функций. Для этого исполь- 
зовалась квадратурная формула, основанная на замене кубическим сплайном слабо осциллирующей части 


2Кт5 
подынтегральной функции, а сильно осциллирующий множитель е *^” рассматривался как весовая функ- 


ция [15]. 
Результаты исследования. Итак, бесконечная пластина лежит на упругом винклеровском основании. 
На нее действует нормальная сила, перемещающаяся по траектории, изображенной на рис 1. 
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Рис. 1. Траектория движения нагрузки 


Расчеты проводились при следующих исходных данных: й = 0,25 м; с=221 м/с; Е= 232469-10 Н/м’; 
у= 0,36. На рис. 2 приведены результаты расчетов, соответствующие К = 1,864 м“ иа= 125м/с. 
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Рис. 2. Изменение перемещений и напряжений 


Данные графики описывают изменения перемещений и напряжений в системе координат, связанной с 
точкой приложения движущейся силы Р. Ось # направлена по касательной к траектории движения {‚ аосьи — 


по внешней нормали к области, ограниченной траекторией п (см. рис. 1) при 2 = р 
При этом вектор перемещений и тензор напряжений соответственно представлялись в виде 
И =Ш-1 п.п +’. , 5=51.И +51. пт+5т-(т+пГ), 


где А — нормаль к пластине. 
На рис. 2 показано изменение вдоль оси { компонент вектора перемещений ИТ, ИИ, Ци! ‚ тензора напря- 
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жений 51, ИЕ, УЕ и изменение этих величин вдоль оси и — Ии,т,Иипи и эт, пи, “ти. 

Согласно расчетам, характер поведения указанных величин во всех точках траектории остается неиз- 
менным. Расчеты показали также, что при изменении скорости движения нагрузки а в пределах от 0 до 125 м/с 
компоненты вектора перемещений и тензора напряжений возрастали незначительно (на 3—4%). 211 
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Для изучения распространения энергии упругих волн вычислялись компоненты вектора Умова — Пой- 
тинга е; =—5, 0; (5, — компоненты тензора напряжений, и; — координаты вектора перемещений), направ- 


ление которого указывает направление распространения энергии, а длина описывает количество энергии, пере- 


носимое через единицу поверхности, перпендикулярной к направлению данного вектора в единицу времени. 
На рис. 3. представлено распространение энергии упругих волн около движущейся сосредоточенной 


силы (положение силы обозначено красной звездочкой). 





Рис. 3. Распространение энергии упругих волн 


Обсуждение и заключения. Анализ полученных результатов показывает, что при указанных выше 
пределах изменения скорости длина вектора Умова — Пойтинга практически пропорциональна скорости движе- 
ния нагрузки и качественная картина распространения энергии вблизи силы слабо изменяется в процессе дви- 
жения. Характер поведения рассчитанных выше перемещений и напряжений во всех точках траектории остает- 
ся неизменным, а их значения слабо зависят от скорости движения нагрузки а при изменении этой скорости в 
пределах от 0 до 125 м/с. 

Применение предложенного метода к бесконечной пластине, лежащей на упругом основании, не ис- 
черпывает его возможности. Он может быть использован при рассмотрении более сложных объектов (плиты 
сложной формы, слоистые плиты, вязкоупругие плиты). Рассмотренный метод отличается от упомянутых выше 
большей экономичностью, так как использует для построения решения уже известные решения задач. Оконча- 
тельное решение выражается в удобном виде — как сумма криволинейных интегралов. 

Полученные результаты могут быть использованы в процессе проектирования дорог. Изучение распро- 
странения энергии упругих волн от движущихся транспортных средств позволит оценить воздействие указан- 
ных волн на строения, расположенные вблизи дороги. С учетом данных о характере изменения перемещений и 
напряжений оценивается износ дорожного покрытия. 
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Введение. В современном производстве при выполнении 
финишных операций важную роль играет центробежно- 
ротационная обработка в среде абразива. Основные пре- 
имущества этого метода отделочно-зачистной обработки: 
высокая производительность, низкая себестоимость и 
широкие технологические возможности. В данном иссле- 
довании рассматривается процесс взаимодействия абра- 
зивной частицы с поверхностью детали в рамках статиче- 
ской контактной задачи теории упругости. При этом учи- 
тывается пластическая деформация в области контакта. 
Материалы и методы. Абразивная частица (корунд) мо- 
делируется линейно упругим телом, модуль Юнга которо- 
го значительно больше, чем у обрабатываемого материала. 
Обрабатываемый материал (сталь) моделируется упруго 
пластическим билинейным телом с применением критерия 
пластичности Мизеса. 

Результаты исследования. Выполнено конечноэлемент- 
ное моделирование рассматриваемых конструкций в САЕ- 
пакете АМЗУБЗ. Смоделирован процесс взаимодействия 
абразивной частицы и поверхности детали, проанализиро- 
вано ее напряженно-деформированное состояние. Пред- 
ставлены результаты численных экспериментов, которые 
позволили установить, как распределяются эквивалентные 
пластические деформации при глубинах внедрения конуса 
0,01 мм и 0,05 мм. Полученные данные, а также области 
значений пластической деформации более | % визуализи- 
рованы в САЕ-пакете АМЗУЪ. 

Обсуждение и заключения. Установлено, что эквивалент- 
ная пластическая деформация пропорциональна глубине 
внедрения (ГВ). Она достигает минимального значения 
0,158 при ГВ = 0,01 мм, максимального 0,825 — при 
ГВ = 0,05 мм. Определены зависимости размеров области 
пластической деформации от ГВ для случаев, когда пла- 
стическая деформация превышает 1 %. При максимальном 
внедрении (0,05 мм) радиус деформации составляет 1 мм, 
глубина — 0,8 мм. На основе данных, полученных в ре- 
зультате проведенного исследования, могут быть выбраны 
параметры технологического процесса (скорость враще- 
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ния, размер абразивной поверхности, масса абразивных 
частиц), которые влияют на взаимодействие между дета- 
лью и абразивной частицей. Рациональный выбор этих 
параметров позволит повысить эффективность обработ ки. 


Ключевые слова: центробежно-ротационной обработки, 
абразивная обработка, контактная задача, пластичность, 
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Введение. В современном производстве при выполнении финишных операций важную роль играет 
центробежно-ротационная обработка (ЦРО) в среде абразива. Основные преимущества этого метода отделочно- 
зачистной обработки: высокая производительность, низкая себестоимость и широкие технологические возмож- 
ности. В данном исследовании рассматривается процесс взаимодействия абразивной частицы с поверхностью 
заготовки в рамках статической задачи теории упругости. При этом учитывается пластическая деформация в 
области контакта. Элемент абразивного зерна в виде усеченного конуса (точнее: круг меньшего диаметра этого 
конуса) взаимодействует с поверхностью детали. При этом следует учитывать трение и пластическую дефор- 
мацию данной поверхности. Кинематические или силовые граничные условия применяются к кругу большего 
диаметра. В случае кинематических условий заданы нормальное и тангенциальное смещения окружности и ее 
вращение. В случае силовых условий задаются сила и момент. Исследуются поля напряжений и эквивалентных 
пластических деформаций вблизи области контакта. 

Моделирование геометрии абразивной частицы центробежно-ротационной обработки. Сущность 
метода ЦРО состоит в том, что абразивные частицы 3 и обрабатываемые детали 4 (рис. 1 [1]) загружаются в 
рабочую камеру и приводятся во вращательное движение вокруг вертикальной оси таким образом, что вся мас- 
са загрузки приобретает форму тора [2, 3], в котором частицы движутся по спиральным траекториям. 
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Рис. 1. Центробежно-ротационная обработка в среде абразива: 
| — неподвижная цилиндрическая обечайка; 2 — ротор; 3 — абразивное зерно; 4 — обрабатываемые детали 


Тороидально-винтовой поток обеспечивается конструкцией рабочей камеры станка, состоящей из не- 
подвижной цилиндрической обечайки | и примыкающего к ней вращающегося дна (ротора) 2, имеющего в рас- 
пространенном варианте форму тарели. Детали 4 загружаются в рабочую камеру навалом вместе с абразивны- 
ми частицами 3. Для уменьшения износа внутренние поверхности дна и неподвижной части рабочей камеры 
покрывают износостойким материалом. Чаще всего используются резиновые или полиуретановые покрытия. 
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Схема построения модели внедрения индентора в виде усеченного конуса предполагает сферическую 
форму абразивной гранулы с набором усеченных конусов [3]. Примерное представление о геометрии сфериче- 
ской абразивной частицы показано на рис. 2 [4]. 


# 





Рис. 2. Геометрия сферической абразивной гранулы: | — приближенная к реальной; 2 — моделируемая 


Согласно модели, абразивное зерно имеет форму усеченного конуса. Такая модель позволяет учесть 
всегда имеющее место затупление абразивных зерен. На рис. 3 представлена схема взаимодействия абразивной 
частицы и детали [5]. 





= 


Рис. 3. Схема внедрения индентора в виде усеченного конуса 
В силе, действующей на индентор, выделяют нормальную и касательную составляющие [2]. Согласно 
работе А. Н. Бескопыльного [5], нормальная составляющая силы резания равна: 
# г* * * 
2ай Е, И 
В 2 
Р= й +С, 
= 


‚ И Й 
С 
* (1 — 2/к) * * 
ле с = ( —_0 } сиё) нь. С =а ( _0 } сиё) ; 0 —_ относительная высота наплыва; 
В 
В приведенный модуль упругости; й” — глубина внедрения зерна; х = 5^от — параметр пластичности. 


Математическая модель. Оптимизация процесса истирания требует разработки усовершенствованных 
моделей фрикционного взаимодействия между абразивными частицами и поверхностью металлического ком- 
понента. В этой модели необходимо учитывать нагрев и износ поверхности из-за удара и скольжения абразив- 
ной частицы. 

При абразивной обработке происходит контактное взаимодействие, приводящее с износу и разогреву 
поверхности обрабатываемой детали. Основные данные о процессах трения и износа представлены в [6]. Обра- 
ботка в ротационной камере описана М. А. Тамаркиным и его учениками [1-5]. Особенностям этого процесса 
посвящены работы [7-9]. 

В представленной статье контактное взаимодействие абразивной частицы и поверхности детали рас- 
сматривается в рамках осесимметричной деформации фрагмента абразивной частицы с поверхностью обраба- 
тываемой детали. Фрагмент абразивной частицы представляет собой усеченный конус, меньшее основание ко- 
торого взаимодействует с поверхностью детали (рис. 4). 
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Рис. 4. Модель частицы и детали: | — абразивная частица; 2 — деталь 


Согласно модели, абразивная частица (корунд) — это линейно упругое тело, модуль Юнга которого 
значительно больше, чем у обрабатываемого материала. Обрабатываемый материал (сталь) моделируется как 
упруго пластическое билинейное тело. При этом применяется критерий пластичности Мизеса [10]. 

Абразивная частица и деталь контактируют без трения. Контактными поверхностями являются верхняя 
плоскость детали, меньшее основание и боковая поверхность конуса. Нижняя плоскость и боковая поверхность 
детали закреплены по нормали. На большем основании усеченного конуса заданы граничные условия: силовые 
(равномерно распределенное давление) или кинематические (вертикальное перемещение). Рассматривается 
внедрение частицы в деталь, исследуются зона пластических деформаций и их максимальные значения. 

Метод решения. Для решения описанных краевых задач используется метод конечных элементов, реа- 
лизованный в САЕ-пакете АМЗУЪ. На рис. 5 представлена конечноэлементная сетка первой задачи. 
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Рис. 5. Конечноэлементное разбиение в контактной задаче: | — абразивная частица; 2 — деталь 


В области контакта и возможной пластической деформации выполнено сгущение сетки. Используются 
конечные элементы РГАМЕ183 (материал 1), РЕАМЕ182 (материал 2) и ТАКОЕ169, СОМТА172 для поверхно- 
стей контакта. Абразивная частица и деталь симметричны, поэтому исследуется половина осевого сечения 


(рис. 5). 
Результаты численных экспериментов. В численных расчетах в задаче использовались следующие 


данные: радиусы усеченного конуса — 0,5 мм и 1 мм; высота конуса — 1 мм; радиус детали — 5 мм; толщина 
— 3 мм. Модуль Юнга материала 1 равен 2х 10° ГПа; коэффициент Пуассона — 0,3. Модуль Юнга материала 2 


равен 2х10° ГПа; коэффициент Пуассона — 0,28. Величина напряжения текучести — 0,22 ГПа. На рис. 6 


представлено распределение эквивалентных пластических деформаций при глубинах внедрения конуса 0,01 мм 
(рис. ба) и 0,05 мм (рис. 60). 
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Рис. 6. Распределение эквивалентных пластических деформаций при глубине внедрения конуса: 
0,01 мм (а); 0,05 мм (6) 


На рис. 7 показаны зависимости максимального значения эквивалентной пластической деформации от 


глубины внедрения (ГВ). 
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Рис. 7. Максимальные пластические деформации в зависимости от глубины внедрения 


В ходе экспериментов по моделированию авторы обнаружили, что эквивалентная пластическая дефор- 
мация пропорциональна глубине внедрения. Эквивалентная пластическая деформация достигает минимального 
значения 0,158 при ГВ = 0,01 мм, максимального 0,825 — при ГВ = 0,05 мм. Определены зависимости размеров 
области пластической деформации от глубины внедрения для случаев, когда пластическая деформация превы- 
шает | % (рис. 8—10). Области значений пластической деформации более 1 % выбраны с помощью программно- 


го обеспечения АМЗУЪ. При этом определены глубина (Н) и радиус (Г.) зоны деформации. 
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218 Рис. 8. Размеры области пластической деформации 
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Рис. 9. Зависимость радиуса области пластической деформации от глубины внедрения 
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Рис. 10. Зависимость глубины области пластической деформации от глубины внедрения 


Определена зависимость глубины внедрения конуса от равномерно распределенного давления, дей- 
ствующего на большее его основание. Эта зависимость представлена в табл. 1. 
Таблица 1 
Зависимость глубины внедрения конуса от давления 


Напряжение (Н/мм”) 200 250 300 





-0,059 —0,119 —0,237 —0,354 -0,530 


Результаты численных экспериментов позволяют определить зоны пластических деформаций И ИХ 


величину в зависимости от глубины внедрения. А данные табл. 1 позволяют связать их с усилием, которое 
необходимо приложить к абразивной частице. Это усилие можно определить по параметрам технологиче- 
ского процесса. 

Обсуждение и заключения. В настоящей работе выполнены: 

— конечноэлементное моделирование процесса взаимодействия абразивной частицы с поверхностью об- 
рабатываемой детали; 
— проведен анализ напряженно-деформированного состояния поверхности. 

Определена зависимость максимальной пластической деформации от глубины внедрения конуса 
(0,ОТ мм - 0,05 мм). Установлено, что эта величина меняется в пределах от 0,158 до 0,825. 

Определены зависимости размеров области пластической деформации от глубины внедрения, для 
случаев, когда пластическая деформация превышает 1 %. При максимальном внедрении (0,05 мм) радиус 
деформации составляет 1 мм, глубина — 0,8 мм. 

На основе этих данных могут быть выбраны параметры технологического процесса (скорость вра- 
щения, размер абразивной поверхности, масса абразивных частиц), которые влияют на взаимодействие 
между деталью и абразивной частицей. Рациональный выбор этих параметров позволит повысить эффек- 
тивность обработки. 
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Введение. Интеграция концепций надежности и 
оптимизации направлена на разработку структур, которые 
должны быть экономичными и надежными. Эта модель 
называется Оптимизация проектирования на основе 
(ВВОО). между 


механическим моделированием, анализом надежности и 


надежности Фактически — связь 
методами оптимизации приводит к очень высоким 
вычислительным затратам и слабой стабильности 
конвергенции. 

Материалы и методы. Для преодоления этих трудностей 
было разработано несколько методов. Методы под 
названием «Индекс показателя надежности» (ВТА) и 
«Метод оценки эффективности» (РМА) - это два 
альтернативных метода. ВТА описывает вероятностное 
ограничение как индекс надежности, в то время как РМА 
был предложен как путь преобразования вероятностной 
меры эффективности. Предложен метод оптимального 
коэффициента безопасности (ОЗЕ) для 
коэффициентов 


расчета 
безопасности, удовлетворяющих 
требуемому уровню надежности, без дополнительных 
вычислительных затрат для 


Уравнения ОЗЕ сформулированы с учетом ВТА и РМА и 


оценки надежности. 
расширены до случая с несколькими отказами. 

Результаты исследования. Несколько линейных и 
нелинейных законов распределения применяются к 
исследованиям композиционных нитей, а затем 
распространяются на множественные режимы отказа. 
Было показано, что идея метода ОЗЕ заключается в том, 
чтобы избежать оценки ограничения надежности с 
помощью конкретного процесса оптимизации. 

Обсуждение и выводы. Упрощенная структура реализации 
стратегии ОЗЕ заключается в разделении анализа 
надежности. Она 


оптимизации И предоставляет 


проектировщикам эффективные решения, которые 


должны экономически удовлетворять требуемый уровень 
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сотраге4 ю ОЗЕ Ваз зеуега! адуааеез: зтаП пишбег оЁ 
орИпитайоп уапа ез, соо сопуегоепсе забИу, зтаП 


сотрийпе Ите, зай$Тасноп оЁ Ше гедие4 тепаб цу 1Теуе|5. 


Кеуумогд$: Кейабщу-Вазед Оез1еп Орнпитайоп, Згасвига| 
КепабИу, Заебу Еас®югз, Еайеие Патазе Апа1у$1$, Ма 
Еаиаге эсепаг10$. 


Еог сйайпоп: КВаптапда Оаз, Ашураз Ппа4 В.. Есет 
орИи заРеёу Расюг арргоасВ Юг зубет тепабу-Базеа 
4ез1еп орниитайоп У аррПсайоп т сотрозйе уагиз. 





надежности. Показано, что КВШО в сравнении с ОЗЕ 
имеет ряд преимуществ: небольшое число переменных 
оптимизации, хорошая стабильность сходимости, малое 
время вычислений, удовлетворение требуемым уровням 
надежности. 


Ключевые слова: оптимизация на основе надежности, 
структурная надежность, факторы безопасности, анализ 
ущерба от изношенности, сценарии множественных 
отказов. 
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1. шогодисйоп 

\Пеп Раегиитзйс Оез1еп Орйплтаноп (ООО) шефо4$ аге изе4, даегииизЯс оритит 4ез1е1$ аге изиаПу 
разВе © Ше 4ез1еп сопзгаше Боцидагу, 1еауше Ш@е ог по тоош Юг юегапсез (ог ипсецашиез) ш 4ез1еп, 
тапа{асвиге, ап орегайп» ргосеззез. Зо даегтилзйс орйтит 4е51еп$ оБМаше4 у/фочЕ сопз14егайойп оЁ 
ипсецашиез соч]4 Теа4 о ипгейае 4ез1ет$, ШегеГоге саШи® Юг Кейабиу-Вазе4 Без1еп Ориплтайоп (КВРО). И 1$ 
Фе оБесйуе о ВВПО ® 4ез12п згасвигез Фаё зВош4 Бе Бо есопопис ап4 тейае. Но\еуег, е соирПи® Бебмееп 
Фе тесвапса| то4деппе, Фе тепаб у апа[узез ап Ше орипитаНоп тефо4$ [еа4$ ю уегу шей сотршанопа| со5 
ап \еаК сопуегоепсе забШиу. То оуегсоште Шезе АЁАсиез, мо рошб оЁ ме\ Вауе Бееп сопз14еге4. Егот а 
гепаб Шу ме\ рошь КВРО шуо|уез е еуаайоп о? ргоба 1 5Ис сопзга, урлсЬ сап Бе ехесие 4 ш хо аШегеп: 
\ауз: ейпег изше Фе Кейаб Ииу ш4ех АрргоасВ (ВТА) ог Фе Регюогтапсе Меазиге АрргоасВ (РМА) [1]. ТВе та]ог 
ЧАсаКу Пез ш Фе еуааноп оЁ Фе ргоБабШ5Ис сопзга, уВлсЬ 1$ ргобПуе[у ехрепууе апа еуеп Чуегое$ Юг 
тапу аррИсаНоп$. Но\еуег, Кот ап орипитайоп ме\ рошь а Ву тефо4 Базе оп зипиЙапеоч$ зоНоп оР ве 
гепаб у ап4 Фе орниитайоп рго ет Ваз зиссеззПу гедисе4 Фе сотршщайопа! ите рго ет. Но\хеуег, Фе вубиа 
ВВОО ргоет 15 шоге сотр]ех ап а ог д&егиитзис 4ез1еп ап тау по Теа4 © [оса| орита. То оуегсоте Бо 
гамфасК$, ап Оритит за у Еасюг (ОЗЕ) шефо4 Ваз Бееп ргорозе4 © сотрие за у Ёасюгз заязРуше а гедлитеа 
гепаб у 1еуе] убой Четап4ше а4Ч1опа| сошрийп® со5 ог Фе тепаб пу еуаайоп. Тре еЁЯслепсу оЁ фе ОЗЕ 
шешо4 Ваз Бееп Четопзиаже4 ге]ануе тю Фе Вуб14 опе Юг а Ппеаг Ч1зитБиноп 1а\м [2]. ш 15 могк, Фе ОЪЕ 1$ 
геГогпи]айе4 сопзАеги>» ВТА апа РМА. М№х ап ежепзоп ® шире №Паге то4ез 1$ сагле оч. ЕтаПу, ап 
аррПсайоп оп сотрозИе уагпз$ 1$ сагле ом юг Ппеаг ап4 попПпеаг а15лБиноп 1а\5 ап шире №Пиаге то4дез 
(зсепа110$). 

2. ВепабШу-Ва$зед Эезеп Орйпи7айоп 
Тре тепабпу-Базе4 4ез1еп орнпитайоп рго ет 15 регогте4 Бу пезИпе мо зиб-ргоЫеитлз [3, 4]: 


1: Орнпихайоп ргоет: ТВе обесйуе 15 ю шшииите ап об]есйуе РшсНоп Х (х) заб] ес{ фо деегиит$ Ис 
опа ©, (х) < 0 ап4д а геаште4 тепаб!иу [еуе] В(х, и) > В, аз ЮПо\5: 

шш: Х(х) 

2: ВепабИиу рго ет: Тве оБесйуе 15 ю Нп4 ше шшипит 41$апсе Бебмееп Фе опош оЁ Фе попта|тед 


зрасе ап4 ве МРР (Моз+ Ргобае Гааге Рош ог 4ез1еп рош®) оп Фе Пий $ае РапсНоп (Ро. 1). ТВе рго ет сап Бе 
\иШеп аз ЮПо\з: 


зиБуесё ю 2) (х) < 0 апа В(х,ч) > В, (1) 


В = па 4 (и) = зиб]есё ю Н(х,а)<0 (2) 


Ш, 
2.4 
1 


Еог тоге 4еа15$ абоиё гейаб у апа[уз$1$, геадегз сап зее [5]. ТВезе мо заб -ргоепл$ аге сагле4 т мо 





Ч1егеп* зрасез: рВуз1са[ ап4 погта|7е4 зрасез (Е1е. 1). 
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Хх Риузтса! 5расе И> № шттайхе4 5 расе 


$ате Еайиге бате Еайиге 


Кезоп 


Кезоп Кезоп 


Дез12п Н(х.и)=0 
Рот А 


и, 





Е!о. 1. РБуз1са| ап погта|7е4 зрасез [3, 4] 
Рис. 1. Физические и нормализованные пространства [3, 4] 


Трезе пефе ргоетз гедшге а меВ сотшрийие име, Во\еуег, ше Оритишт За у Еасюг (ОЗЕ) 1$ изе4 ф№ю зо|уе 
Фезе питегса1 АТАсш вез. 

3. ОЗЕ деу@ортеп 5 

п Фе ВТА апа РМА, Фе тауог АЁАси Гу Пез ш е еуаайоп оЁ Фе ргобаб зас (тейабИиу) сопзгаи 5; 15 
1$ рго1ауеу ехрепзуе ап4 еуеп Фуегое$ Юг тапу аррПсаНоп$ [6]. ВТА 4езстез$ @е ргоБаБ1\5ис сопзгатиЕ аз а 
гепаб у шдех “бе РМА \аз ргорозеа Бу сопуегте Фе ргобаб у теазиге 0 а ре!югтапсе теазите. [ Ваз Бееп 
Кип ш Фе Шегабиге фа РМА 15$ тоге еРйслет{ апа а е ап КТА ш Фе КВПО ргосез$ [1]. Еифегтоге, РМА 
сопуегоез ю ЧШегеп оритит зо[айоп$ упеп загиие Нот @Шегепе п1йа1 4ез1епз [7]. п 1$ зесйоп, и 15 зпо\’и фа% 
фе ОЗЕ юопищаНоп$ сап Бе дедицеа гота Бо арргоасВез$ (ВТА & РМА). 

3 ОЗЕ сопу4Аегто ВТА апд РМА 

3.1. ОЗЕРЕ разе4 оп ВТА 

Тре ргобаб!1$с сопзмашЕ ш ВТА сап Бе еуашаеа4 Ъу зо]уше Фе Игз{-ог4аег гейа Шу апа[у$1$ ю сасяае 
Фе 4ез1еп рошё Р*, уме 1$ Югишае4 аз ап ориптайоп рго ет у ап едиаШу сопзгашй Н (и) = 0, ш Фе 
погта|7е4 расе: 


шш : (а) 
т (За) 
зиб]есёю: Н(и)=0 
апа сап Бе а[50 \и@еп аз: 
шш : а“ (и) 
и. (ЗБ) 
зиб]есёю: Н(и)=0 
Тре Гастапелап Раисйоп Юг Фе рго ет (3Ъ) сап Бе уп@еп а$ 
Г(и,^.) = а* (а)+^-Н (а) (4) 
Тре оритаШу сопд1оп$ ог Фе Гастапелап Гапсйоп аге: 
2 
ОЕ _ 04 `. ОН _ | ‚ 1=1,...П (>а) 
ди; и; ди; 
ОГ, 
=. =Н (в) (55) 
Озше \е ехргезз10п Юг 1е запаге Ч апсе 4 ш 1е едиайоп (35), ме се: 
Н 
Е: п АН (6) 
2 ди; 


3.1.2 ОЗЕ Базе оп РМА 
Тре еуааноп ог Ше ргобаб1$с сопугаш ш РМА гедиез$ ап шуегзе гейаб (у апа|уз1$, УмсЬ соггезроп@$ 
{0 Фе шуегзе рго ет оЁ Фе гейа Шу апа1уз1$ 10 саси=жще 1е 4ез1еп рош Р*, \шсь 1$ огиа[ае4 аз ап орниитайоп 


ргоет \И® ап едиаШу сопзаате а? (и) = В: ‚ ш погта|7е4 зрасе: 
и] : Ноа) 


Та 
зибуеси ю: 4 (а) =В, р 
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ап4 сап Бе а[$0 ууп@еп аз: 


шш ‚: На) 
ь (76) 
забуесе 0: 4“ (и) = В 
Тре Гастапелап Раисйоп Юг Фе рго ет (76) сап Бе уп@еп а$ 
Е(а,^,5) = Н(а)+^. 4? (и) В | (8) 
Тре оритаШу сопд1оп$ г Фе Гастапелап Гапсйоп аге: 
2 
Е _ ОН ›, 94 _| - 1=1[..... И (Эа) 
ди; ди; ди; 
ОГ 2 
—=а (и) =0 9Б 
О), } ты 
Озше \е ехргезз10п Кг 1е запаге Ч апсе 4 ш 1е едцайоп (76), \уе се: 
и; Но: И (10) 
2), ди; 


Ассог4ше 1 еацаНоп$ 6 ап 10, Ше тези@те погта|те уесюг и сап Бе уийеп ш {ети оЁ Фе зате 
дегуануе ОН /ди, . 

3.2 ОЗЕ Гог сотропет ВВОО 

[её и5 соп$14ег по\’ Фе сазе оР и=2? попта|те уапма Без г = 1, 2 (5ее Е1е. 2). ТВе фапееп оРо1$ элуеп 
Бу: (апа = и, / м. Озшх (6) апа (10), ууе се Фе зате Гогтшайоп$ аз юПо\з: 
ОН 
био 
ОН 
би 


Ир 


ап а = 


Тре Юги\щайоп (11) зВо\/з фай Феге 15 по АШегепсе Бебмееп ВТА апа РМА Юг Фе ЮПомлие ОЗЕ деуеюртеп5. 
Тре ргоеплаз$ (3) апа (7) уе ч$ Фе гепабИиу шаех В аз Ме тшипит 41${апсе Бебмееп Фе Пти зае зиГГасе апа фе 


опр [8]. ТЬ1$ шеап$ фай Фе гез пе гепаб у ш4ех тау Бе 1о\уег ог Б1еБег ап Ше {агое( гепаб Шу шдехВ,. Аз 
\е ул {о за1$Ру а гедииге4 фагоей гепаБ у 1еуе] Гог Фе орипл1тайНоп рго ет, \уе сап \тие 
И 
2 2 
=» и (12) 


1=1 


Таз, Бе сепега] ехргеззлоп Гог Ше погта|7е4 уапа е и; уВеп изшо ВТА ог РМА, 15 


2 
он 
[2 


5") 


Л 





2=1,...п (13) 





Тре са1сфаноп оЁ Фе погта|7е4 отаФепЕ ОН / ди 1$ по{ 4иесЙу ассезз1]е Бесаизе Фе тесватса| апа[уз1$ 1$ 
сагле4 ой ш Фе рБуз1са| зрасе гаег ап ш Фе запдага зрасе [9, 10]. Но\жеуег, изше Феогу оЁ %айз$Ис$, \е сап 
депуе фе юоПоулие ехргезулоп Нот \мсЬ Фе сопршайоп оЁ Фе погта|7е4 огаФеп{ сап Бе сагле4 оц Бу арр|у® 


Бе сваш ге оп Фе рБуз1са| огаЧет ОС / ду: 


дн _ дб 07; (х,и) 


Ари. РР (14) 
би О; ОИ 


У\Пеге Т’(х,а) 18 Фе шуеме@ ргобабИзис напзЮгтаной ЕапсНоп. \Ме Ни Фа е поглаН2е4 отафепе сап Бе 


ехргеззе4 аз 
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—- ы ый (15) 
ОИ, ду; 


То зитанете, \е зайзРу Фе гедлиге4 {агое{ гепаб Шу шдех аз №ЮПо\з: 





‚=1,...й (16) 





ТБаз, ме сошрще Фе оритит заеу Расюг$ заб]есё ю Фе оритишт уа[чез оЁ Фе погта|7е уампа ез и; . Ассог4ше 


© (13) апа (15), И Ваз Бееп дЧетопзгае4 Ша Фе оритит уаез оЁ и; сап Бе утШеп ш фе ЮПо\лие Юг: 





р а В. (17) 








уУ*Пеге е 51оп оЁ= 4ерепа$ оп Ве $1еп ое 4епуайуе, 1.е.: 


О беы, 50: 00 Зои 20 „ая (18) 


1 1 





Е1с. 2. Оез1еп рот тодейпе ш а доче погтла|7е4 уапаЫе зрасе 
Рис. 2. Моделирование точек проектирования в двойном нормализованном переменном пространстве 


Тре еЁйслепсу оЁ Фе ОЗЕ 15 © Нап$Рег Фе тепабИиу рго ет шю ап апа[уйса| огищайоп (17). ТЬ$ 
Гоги] айоп Ваз Бееп ежепае4 тю зеуега| 41зитбийоп [а\уз (погиаа1, [оепогипа1, ип огт Уеой, апа Сите] 1а\ $) Юг 
ше аПаге то4е сазез сопз14егие оту ВТА [11]. Ножеуег, ш 1$ \огк, Бо ВТА апд РМА аге изе4 {ю зВо\и Ше 
еЁИслепсу ог Фе ОЗЕ Югтчайоп, ап4 ап ежепзлоп 10 шире ЁаИте сазе 1$ саге ом ш Фе пехё зесйоп. 

3.3 ОЪЕ Гог зует ВВОО 

Тре зубет гепабшу ргоет сап Бе уи6еп аз: 


: _ 2 2 
Взужет = пи Ч(и;) = Дир +и> +...+и, 54: Н, (чмо ,...ии ) < 0 (19) 
уПеге (и; ) 15 Фе шшипиат 41${апсе Бебмееп фе Моз( РгобаЫе Гайоге Ролш( (МРР) ап4 Фе орита| зошоп 
ш Фе поппа|2е4 расе (Е12. 1). АпаН, (и,) <0 гергезеп Фе Чегет Р№Пиаге тодез. Озше зипПаг \ау, Фе 


апа|уйса| огти]аноп изше ОЗЕ Юга зует КВШО сап Ъе уеп аз оПо\з: 
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(20) 








Тье огту|айоп (20) ргоу1Аез ЧШНегет{ оритит уа[аез оЁ Ве погта|7еа уапа ез аё Фе ТаПиге рош{ апа 


{акте шо ассоипЕ зеуега1 Гааге то4е$. 

4. Митетса1 Арр|сайоп оп сотрозе уагп$ 

Сотроз йе таёега15 Вауе Бееп изе4 ш збгасфагез Юг сепиме$. ш гесепе Итез, сотрозЦе раз Вауе Бееп 
изе4 ежепзлуе[у ш аисгай згасвигез, аюто Без, зрогиие 2004$, ап4 тапу сопзитег ргодис&5. Сотрозйе табета[5 
аге Бозе сомаште тоге ап опе Бопае4 шаега], еасв ул а 1Шегепй эбгистага| ргорегиез. Тве таш адуащазе оЁ 
сотрозе шаета[$ 15 фе роепвйа[ ог а 612 гайо ог зи Нпез$ ю \е12Е. То зо]уе Фе сотрозЦе уагиз згисагез (ог 


5у3{еп15), а $тез$ ип Иаега| 5би4у 15 гедлигед. 





а) Корез оРсотрозКе уагп$ Ь) зесноп ое за е4 горе 
Е1е. 3. Мжепа| згасаге 


Рис. 3. Структура материала 


4.1 Ргоет дезсириоп 
Еог и1-таёема| згасбиге Шизгае ш Е1э. 3, фе патбег оЁ уаги$ п, , Ше зесноп 5,;, фе Уоцпз'5 шоди$ Е), 
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Рис. 4. Напряжение действующие на нити пряжи 
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4.3 Ргоседиге$ 
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КВРО ргосейиге 
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гепабу-базе4 оритит 4е51еп$ Уфо ад 11опа! сотрийпо со$ё Юг ргоба15ас (тепаб у) сопзгаш ап4 сап 
|еа4 №ю 21оба! орита. [ 15 зВо\и Фа Фе КВРО сошрагед ю ОЗЕ Баз зеуега|] адуатасез: зтаП пипабег оЁ 
орНиитайоп уапа е$, соо сопуегоепсе за Шу, эта сотрийпе ите, зай$РасНоп оЁ Фе гедиеа тепабиу 1еуе|5. 
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Введение. В настоящее время на производствах с объемом 
выпуска до 4000 булок в день используются 
тестоделительные машины вакуумного типа. В процессе 
эксплуатации тестоделительного устройства из-за износов 
рабочих поверхностей поршня, камеры, барабана зазор 
между ними выходит за величину, равную 50 мкм, при 
которой обеспечивается вакуум во всасывающей камере. 
В результате этого процесс всасывания становится 
нестабильным, тестоделительное устройство нарушает 
точность развесовки хлебобулочных изделий. Ремонт 
такого оборудования проводится в основном с 


использованием полной или частичной замены 
изношенных деталей и узлов на новые. Для повышения их 
долговечности возникает потребность в разработке новой 
высокоэффективной технологии с восстановлением 
изношенных поверхностей деталей сварочно- 
наплавочными методами. 

Материалы и методы. Представлена новая методика 
определения количества объектов для исследования с 
использованием программы «Статистика». Определены 
поверхности износа деталей тестоделительных машин 
вакуумного типа. 

Результаты исследования. Проведены микрометражные 
исследования деталей  тестоделительных устройств, 
которые показали наличие у них значительных искажений 
размеров из-за локального износа рабочих поверхностей. 
При этом боковой зазор между всасывающей камерой и 
главным поршнем, барабаном и всасывающей камерой в 6 
раз превышает допустимый, а вертикальный зазор между 
мерной камерой и поршнем превышает допустимый зазор 
более чем в 10 раз. Износы рабочих поверхностей деталей 
тестоделительных машин носят локальный характер, при 
этом диапазон значений составляет: для главного поршня 
— 10-200 мкм; мерного поршня — 250-900 мкм; 
всасывающих и мерных камер — 300-400 мкм; 


поверхностей барабана — 280—300 мкм. 
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Обсуждение и заключения. Проведенные Оксиубот апа Сопсибоп5. ТВе сопдище писготейлс 


И ОХ УЧЕНЬЕ ЕЩЕ ЗЛА ПОЗА НТ заЧ1ез зпо\е4 Ше ргезепсе оРарргесла е °12е Ч1зюгНоп$ 4ие 


стройств показали наличие них значительных 
устр у {о Ше 1оса| \еаг оЁ Фе уогкше загсез. Вазе4 оп Фе гези 5 


искажений размеров из-за локального износа рабочих 


ы Ш 1 п Г |] 1 
поверхностей. Основываясь на полученных результатах, оаште4, 1 сап Бе агвие4 1паё Ше тозё ргофиснуе ап 


можно утверждать, что наиболее технологичным и есопописаПу улае 1есбиаие Юг Фе тезогайоп оЁ \от 


экономически целесообразным для восстановления асе оЁ доче аглаег раг5 15, ЮГ ехатр[е, {фе 


изношенных Поверхностеи деталеи тестоделительных е]есвгозрагК тас1и1по. 


машин является метод электроискровой обработки. 


Ключевые слова: тестоделительное стройство, Кеууог4$: доиор 4у14ег, дочоВ ууелоБ{ 4151иНоп, \еаг оЁ 
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развесовка теста, износы деталей тестоделительного Дои2В Д1у14ег раз, теазигетейе Фастат, гарз$ ш тобШе 
устройства, схема измерения, зазоры в подвижных ппегГасез оР дои2В ату1дег. 


сопряжениях тестоделительных устройств. 
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Введение. Хлебопечение в сфере пищевой промышленности занимает особое место. Продукция 
предприятий этой специализации является наиболее востребованной. Однако, несмотря на большой спрос на 
хлеб и хлебобулочные изделия, подобные производства в некоторых случаях могут быть убыточными. Связано 
это с недостаточной механизацией процесса выпечки и трудоемкостью их обслуживания. Чтобы сделать 
предприятие подобной специализации рентабельным, необходимо использовать в процессе приготовления 
хлеба высокотехнологичное оборудование, к которому относятся, например, тестоделительные машины [1]. 

Тестоделительные машины предназначены для механической обработки теста путем деления его на 
порционные куски определенного веса и формы. Используются тестоделительные машины, чаще всего, в 
пекарнях, поставляющих на рынок хлеб, батоны, булочки ит. д. 

Наиболее широкое распространение на хлебозаводах и пекарнях получили тестоделительные машины 
отечественного производства «Восход», «СОЮ/ЛОМТ» и импортного — «РАКТА Ц2», «Киткауа», «СВТ» [2, 3]. 
Импортное оборудование, хотя и является более функциональным, чем аналоговые отечественные модели, 
однако стоит в 1[,5—2,5 раза дороже. 

Главным узлом тестоделительной машины, выполняющим основную функцию, является делительное 
устройство (рис. 1), которое состоит из бункера А; воронки бункера В; корпуса всасывающей камеры С; 
всасывающей камеры КЁ; делительного ножа 0); главного поршня ЕЁ; мерного поршня С; задней заслонки / 
барабанного типа (в дальнейшем — барабан) с мерной камерой Н [4]. 
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Рис. 1. Устройство тестоделительной машины вакуумного типа: А — бункер; 
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В — воронка бункера; С — корпус всасывающей камеры; О — делительный нож; Е — главный поршень; 
Е — всасывающая камера; @ — мерный поршень; Н — мерная камера; [Г — задняя заслонка (барабан); 
К — конвейерная лента 
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В процессе эксплуатации изнашиваются рабочие поверхности главного поршня (рис. 2 а), корпуса 
всасывающей камеры (рис. 2 6), барабана (рис. 2 в) и др. 


А 





Рис. 2. Общий вид деталей тестоделительного устройства 
(А — зона максимального износа): а — мерный поршень; 
© — корпус всасывающей камеры; в — барабан 


Вследствие этого зазоры между ними выходят за величину, при которой обеспечивается вакуум во 
всасывающей камере, процесс всасывания становится нестабильным и тестоделительное устройство нарушает 
точность развесовки хлебобулочных изделий [5]. 

Ремонт такого оборудования до настоящего времени на 70-80 % проводится с использованием 
запасных частей стоимостью от 200 до 300 тыс. руб. Особенно высоки цены на комплектующие детали и узлы 
оборудования зарубежного производства [6]. 

На ряде предприятий технического сервиса восстановление деталей проводится методом наращивания 
чугуна, нержавеющей стали и цветного металла на изношенные поверхности, а затем механической обработки, 
согласно допусков и технических требований завода-изготовителя [6, 7]. Однако данный метод экологически 
небезопасный, трудоемкий и имеет высокую себестоимость. 

Для повышения долговечности таких машин необходима новая высокоэффективная технология, 
позволяющая восстанавливать детали нанесением покрытий с заданными физико-механическими свойствами. 
Одним из основных критериев выбора метода восстановления изношенных деталей является оценка значений 
износов поверхностей деталей, определяющих необходимую толщину металлопокрытий. 

В связи с этим целью работы является оценка средних значений зазоров и износов деталей 
тестоделительных устройств по результатам микрометражных исследований. 

Материалы и методы. Выбор количества исследуемых тестоделительных устройств для проведения 
микрометражных исследований проведен с использованием критерия хи-квадрат. При использовании данного 
критерия задавались критическим значением мощности р» = 0,80 [8] и значением односторонней 


доверительной вероятности р„ = 0,80 [9]. 


Определение количества объектов производится по выдвинутой гипотезе: при заданном значении 
уровня значимости & = 0,05. По нулевой гипотезе количество объектов достаточно для исследования, по 
альтернативной — нет. Если уровень значимости текущего значения @а; выше принятого значения 0,05 и 
текущее значение мощности критерия хи-квадрат р»; выше критического значения 0,80, то имеет место нулевая 
гипотеза, и наоборот. 

Для принятых значений параметров график зависимости количества объектов № от односторонней 


доверительной вероятности р„ критерия хи-квадрат представлен на рис. 3. 
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Рис. 3. График зависимости количества объектов М от односторонней доверительной вероятности 


р. критерия хи-квадрат 


Из графика видно, что для принятой односторонней доверительной вероятности р, = 0,80 количество 


объектов для проведения эксперимента равно М№ = 19 шт. 

Результаты расчетов в программе «5ай5Иса» показывают, что расчетное значение уровня значимости 
критерия хи-квадрат а; = 0,064 выше принятого 0,05, а действительная мощность критерия р,; = 0,84 выше 
критического значения 0,80 [8]. Полученные результаты подтверждают нулевую гипотезу при определенном 
количестве тестоделительных устройств М = 19 шт. 

При проведении микрометражных исследований на первом этапе замерялись зазоры в сопряжениях. 

Согласно рис. 1, тестоделительное устройство состоит из всасывающей и мерительной камеры, внутри 
которой перемещаются главный поршень и мерительный соответственно. 

При измерении зазора в горизонтальной плоскости этих соединений поршень | сдвигается до упора к 
неподвижной камере 2 и в противоположной стороне измеряется зазор // в разных сечениях (рис. 4 а). Также в 
этом положении поршня 1 измеряется зазор 12 в вертикальной плоскости. 





Рис. 4. Схема измерения зазора в подвижных соединениях 


деталей тестоделительных устройств 


Затем поршень | сдвигается в противоположную сторону (рис. 4 6) и алгоритм измерения повторяется. 
Такие измерения проводились в крайних и среднем рабочем положении главного и мерительного поршня 1 
(рис. 1). 

Соединение задняя заслонка-барабан конструктивно выполнено из двух частей — корпуса с мерной 
камерой и накладки. Зазор в данном соединении измеряется между барабаном и всасывающей камерой. 
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Если зазор в соединениях больше 0,05 мм, измеряется износ рабочих поверхностей деталей. При работе 
тестоделительного устройства в соединениях главный поршень — всасывающая камера, мерительный поршень 
— мерительная камера и барабан — всасывающая камера (рис. 1). Неизношенная поверхность деталей 
отсутствует, их относительный износ измеряется согласно [10, 11]. 

Относительный износ наружных поверхностей деталей в 1-м сечении вычисляется по формуле 

Н 
и: № Чтлах — Ч, 
где 4 и.х — максимальный размер детали, мкм; 4, — размер детали в 1-м сечении, мкм. 
Относительный износ внутренних поверхностей деталей в 1-м сечении вычисляется по формуле 


и ыы р; —- а 


где 0); — размер детали в 1-м сечении, мкм; Оии — минимальный размер детали, мкм. 

Сечения деталей камер тестоделительных машин при измерении выбираются согласно рабочему ходу 
поршней в тех местах, где контактируют рабочие поверхности (рис. 1). 

Относительный износ поверхности барабана измеряется путем радиального отклонения 
фиксированных точек поверхности барабана относительно центровых отверстий цапф подшипниковых опор. 


Схема измерений рабочей поверхности барабана приведена на рис. 5. 
1 2 3 4 56 7 89 
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Рис. 5. Схема измерений рабочей поверхности барабана 


Для измерений барабан устанавливается в центрах токарного станка. Измерения при продольном 
перемещении производятся индикаторной головкой вдоль барабана в разных плоскостях, но без изменения 
положения этой головки в радиальном направлении. 

Для измерения наружных деталей тестоделительных устройств используется микрометр МРО2102 с 
погрешностью 0,003 мм, ценой деления 0,002. Для внутренних деталей — индикаторный нутромер НИ-100М 
ГОСТ 868, погрешность измерений которого составляет 0,015 мм, с индикаторной головкой часового типа 
марки ИЧ1О ГОСТ 577, класс точности 1, с ценой деления 0,01 мм. 

Результаты исследования. Результаты средних значений зазоров между поршнями и камерами 
тестоделительного устройства, согласно рис. 3, приведены в таблицах 1 и 2. 

Таблица 1 
Средние значения зазора между всасывающей камерой и главным поршнем 
Расстояние от переднего Среднее значение зазора, мм (+0,02) 
края камеры, мм 
> > 
> > 
> > 


1 
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Таблица 2 


Средние значения зазора между мерной камерой и поршнем 


Расстояние от переднего Среднее значение зазора, мм (+0,02) 


края камеры, мм 7 
3 





Из таблиц | и 2 видно, что боковой зазор между изношенными всасывающей камерой и главным 
поршнем достигает на отдельных участках значения 300 мкм. С учетом повышенного износа выходной части 
нижней плоскости камеры вертикальный зазор составляет не менее 400—500 мкм. Значения бокового и 
вертикального зазоров между изношенными мерными камерами и поршнями также значительно превышают 
допустимую величину и достигают 550 мкм. 


Средние значения зазора между барабаном и всасывающей камерой приведены в таблице 3. 


Таблица 3 
Средние значения зазора между барабаном и всасывающей камерой 
НЯ Значения зазора в мм (+0,02) в точках 
ИИ Е ЗЕ ВЕ 


в | <05 | 026 | 030 | 050 | 058 | 0.25 | 025 | 028 <05 


в | <0.05 [| 047 [| 015 | 08 | 0 вв | 05 | 05 [055 


до о ож [о [0 [05 [0% | 008 [4005 





Как следует из таблицы 3, величина зазора между барабаном и всасывающей камерой в 300 мкм 


значительно превышает допустимую величину. 
В таблицах 4—5 приведены схемы измерения и средние значения износа поверхностей деталей, 


определяющие работоспособность тестоделительных устройств. 


Мартынова Е. Г. и др. Результаты микрометражных исследований деталей тестоделительных машин вакуумного типа 
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Таблица 4 
Результаты микрометражных измерений главного поршня 


Средние значения 
Сечения .. Наблюдения 
измерений, мм 
3 


Измерения высоты поршня по сечениям, +0,01 мм 








Средние значения 


измерений, мм 





99,95 В зонах 1-2 на нижней поверхности износ глубиной до 0,2 мм на 
99,93 длине 20 мм шириной 5—6 мм 
99,90 


99,86 

99,96 
задиры и царапины глубиной до 0,1 мм, 
10 99.87 отдельные вмятины 
е 


13 99,84 


14 99,80 


Измерения ширины поршня по сечениям, +0,05 мм 


На всей площади нижней и верхней поверхностей зон 3-14 





— 
— 





Средние значения измерений, мм 





2 189,65 
| 
На всей площади боковых поверхностей задиры и царапины 
глубиной до 0,1 мм 
р 19 
р 19° 
р 19 


10 189,80 В этой зоне на правой стороне продольные задиры 
глубиной до 0,3 мм 
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Таблица 5 
Результаты микрометражных измерений мерных поршней 


Результаты измерений, мм 
Сечения Наблюдения 
(-0,01 м 


Средние 1 2 3 Схема 


значения измерений 


10 


измерений, мм ширины 


11 


поршня по 
12 


—— 
=— 
ОИ 


сечениям 


ПОРШЕНЬ №1 (левый), измерения высоты поршня по сечениям 


99,88 
На нижней поверхности в зонах 1-3 задиры 
99,88 и 
глубиной до 0,3 мм 
На верхней поверхности (половине площади) в зонах 
99,40...50 10-12 усталостный износ, на остальной — задиры 

99,55 глубиной до 0,2 мм 

ПОРШЕНЬ №1 (левый), измерения ширины поршня по сечениям 


82,33 
82,33 
О. 
82,38 
82,36 
82,31 


81,92 
82,00 
82,15 


81,80...90 
82,00 


82,14 








В зонах 7-9 однородный износ, величина износа левой 
поверхности больше, чем правой 


. 2 
В зоне 10 усталостный износ площадью до 3,5 см‘. 


Однородный износ левой поверхности больше, чем 
правой 
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ПОРШЕНЬ №2 (правый), измерения высоты поршня по сечениям 


рт 

О 99, На нижней поверхности в зонах |-6 задиры и царапины 
ОЕ ПО - ОИ 

О О УСО 

О ОО ООО ООО 
о ПОНИ 
На верхней поверхности в зонах [0-12 усталостный 
износ на площади до 90% 


99,20...50 


ПОРШЕНЬ №2 (правый), измерения ширины поршня по сечениям 


82,29 
82,31 
82,18 
82,38 
82,41 
82,28 
81,66 
81,98 
82,22 


81,49...75 На обеих поверхностях в зонах 10-11 усталостный 


$1 58 9 износ на площади до 50% 


Как следует из таблицы 4, износ высоты главного поршня по сечениям неравномерен. Он минимален в 


дл > 


— 


задней части поршня (10-20 мкм) и постепенно увеличивается до максимального значения 150 мкм к его 
переднему торцу. Аналогичен характер изменения размеров поршня по ширине, у которого максимальное 
значение износа не менее 200 мкм также находится у его переднего торца. На всей площади поверхностей 
трения видны борозды различной глубины (от 100 до 300 мкм), царапины, отдельные вмятины. Износ поршня 
связан с механическим истиранием. 

Результаты обследования мерных поршней (таблица 5) показывают, что износ по высоте и ширине 
максимален в передней его части (таблица 3) и может достигать 800—900 мкм. Однако имеется особенность, 
вызванная наличием вертикального и горизонтального окон. Участки трущихся поверхностей, находящиеся 
против этих окон, изношены больше — в среднем на 250-300 мкм. Кроме того, в передней заходной части 
поршней, напротив окон, на верхней и боковых поверхностях имеются участки повышенного износа. Это 
отчетливо видно на рис. 2 а. 

Исследование износов деталей всасывающей и мерной камер показало, что длительная эксплуатация 
также отразилась на изменении их размеров. Так же, как и поршни, они изношены неравномерно по сечениям 
по всей длине рабочего хода. Характерным для всасывающей камеры является повышенный износ (300-— 
400 мкм) нижней плоскости у выходной части по всем сечениям (рис. 2 6). Ширина полосы повышенного 
износа всасывающей камеры составляет 25—30 мм. 

Результаты измерений радиального отклонения рабочей поверхности барабана, выполненные согласно 
схеме рис. 4, приведены в таблице 6. 
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Таблица 6 


ы | 
Средние значения радиального отклонения рабочей поверхности барабана 


Результаты измерений по сечениям 
Сечения 


Примечание. 


1) >> и < — направления продольного перемещения индикаторной головки; 
2) точки Ж1 и Ж9 - базовые. 





Наряду с измерением радиального отклонения фиксированных точек поверхности, можно, пользуясь 
данными таблицы 3, определить границы изношенных участков и значения износа. Наиболее изношенные 
участки рабочей поверхности барабана находятся в сечениях напротив мерных камер в верхней и нижней части 
накладки (рис. 2 в, 4). 

Как следует из таблицы 3, максимальный износ на наиболее изношенных участках рабочей 
поверхности барабана достигает 280-300 мкм. 

Для восстановления изношенных поверхностей деталей тестоделительных машин в условиях 
современного ремонтного производства, наиболее технологичным и экономически целесообразным является 
метод электроискровой обработки [12]. 

Электроискровая обработка — это процесс перенесения материала электрода на обрабатываемую 
поверхность искровым электрическим разрядом. Особенностью данного метода является локальность 
обработки изношенных участков, отсутствие нагрева деталей и возможность нанесения любых 
токопроводящих материалов, обеспечивающих высокие триботехнические свойства пар трения, при 
минимальной подготовке поверхностей. 

Заключение. Таким образом, проведенные микрометражные исследования деталей тестоделительных 
устройств показали наличие у них значительных искажений размеров из-за локального износа рабочих 
поверхностей. 

Исследования зазоров в соединениях тестоделительных машин показали: боковой зазор между 
всасывающей камерой и главным поршнем в 6 раз превышает допустимый и достигает при среднем значении 
300 мкм; боковой и вертикальный зазоры между мерной камерой и поршнем составляют 550 мкм, что более 
чем в 10 раз превышает допустимый зазор; зазор между барабаном и всасывающей камерой в 6 раз превышает 
допустимый зазор и достигает при среднем значении 300 мкм. 

Износы рабочих поверхностей деталей тестоделительных машин носят локальный характер. При этом 
диапазон значений составляет: для главного поршня — 10-200 мкм; мерного поршня — 250-900 мкм; 
всасывающих и мерных камер — 300—400 мкм; поверхностей барабана — 280-300 мкм. 

Основываясь на полученных результатах, можно утверждать, что наиболее технологичным и 
экономически целесообразным для восстановления изношенных поверхностей деталей тестоделительных 
машин является метод электроискровой обработки. 
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Теоретические основы расчета системы управления гидравлического привода стенда 


ДлЯ испытаний поршневых гидравлических цилиндров” 
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Введение. Долговечность и работоспособность гидравли- 
ческих машин определяется в результате ресурсных испы- 
таний. При этом для силового нагружения гидравлическо- 
го двигателя применяются различные тормозные устрой- 
ства (механические, электрические, гидравлические и др.), 
в результате чего теряется значительное количество энер- 
гии. Этого можно избежать, если при ресурсных испыта- 
ниях использовать метод вращательного движения с реку- 
перацией энергии. Такой подход применим для гидравли- 
ческих насосов, моторов, а также гидравлических цилин- 
дров. 

Материалы и методы. Представлен испытательный 
стенд, конструкция которого позволяет воссоздать усло- 
вия, максимально соответствующие реальной эксплуата- 
ции гидравлических цилиндров. При этом возможна реку- 
перация энергии. Для решения задач исследования ис- 
пользованы методы математического моделирования, 
рассчитаны основные функциональные параметры пред- 
лагаемой конструкции. Определение приращения давле- 
ния в различных точках гидравлической системы базиру- 
ется на теории объемной жесткости. При моделировании 
движения подвижных элементов гидравлической системы 
стенда использованы законы движения ротора. 
Результаты исследования. В структуре испытательного 
стенда рассматриваемые гидроцилиндры размещены в 
напорной магистрали между гидронасосом и гидромото- 
ром. Это позволяет существенно уменьшить сам стенд и 
сэкономить значительное количество энергии за счет ее 
рекуперации. Приведена принципиальная гидравлическая 
схема стенда для испытаний поршневых гидроцилиндров, 
в рамках которой показана работа подвижных элементов 
системы. Выполнено математическое моделирование гид- 
равлической системы стенда. Показана кинематическая 
схема механизма передачи движения между испытуемыми 
цилиндрами. 

Обсуждение и заключения. Представленная в статье си- 
стема уравнений показывает, каким образом определяется 
приращение давления в выбранных узловых точках си- 
стемы рекуперации энергии (в частности, как приращение 
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зависит от времени, приведенного коэффициента объем- 
ной жесткости, расхода рабочей жидкости, площади 
поршней). Величины скоростей перемещения плунжеров 
гидравлических цилиндров определены согласно кинема- 
тической схеме механической передачи стенда. Итоги 
исследования позволяют утверждать, что, благодаря пред- 
ставленному в статье решению, результаты ресурсных 
испытаний гидравлических цилиндров будут адекватно 
отражать их работу при номи нальном режиме эксплуата- 
ЦИИ. 


Ключевые слова: поршневые гидравлические цилиндры, 
стенд, испытания, рекуперация энергии, математическое 
моделирование, кинематическая схема передачи движе- 
НИЯ. 
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Введение. Одним из важных этапов машиностроительного производства, в том числе и производства 
гидравлических машин, является испытание конечного продукта на соответствие техническим требованиям [1]. 
трудо- И 
—щ ресурсные. Они позволяют определить долговечность, работоспособность оборудования в процессе дли- 


Наиболее значимые (при этом энергозатратные) испытания 
тельного применения по назначению. Ресурсные испытания следует проводить в режиме, максимально при- 
ближенном к номинальному режиму эксплуатации гидравлических цилиндров. При этом для силового нагру- 
жения гидравлического двигателя применяются различные тормозные устройства (механические, электриче- 
ские, гидравлические и другие), в результате чего теряется значительное количество энергии, переходящей в 
тепло. Особенно это касается испытаний гидравлических машин средней и высокой мощности. 

В результате активных поисков решения обозначенной проблемы разработан метод ресурсных испыта- 
ний гидравлических машин вращательного движения с рекуперацией энергии [2—4]. Такой подход дает значи- 
тельную экономию при испытаниях гидравлических насосов и гидравлических моторов. Для гидравлических 
цилиндров также разработаны способы испытаний с рекуперацией энергии [5—10]. Схемы, описанные в [9, 10], 
позволяют изготовить стенд, обеспечивающий полное соответствие режима функционирования гидроцилин- 
дров их работе в реальных условиях эксплуатации. 


Материалы и методы 


Постановка задачи. Предложена конструкция стенда, которая позволяет в ходе испытаний воссоздать 
условия, максимально соответствующие реальной эксплуатации гидроцилиндров. При этом возможна рекупе- 
рация энергии, что значительно снижает ее затраты, особенно при ресурсных испытаниях. Для решения задач 
исследования использованы методы математического моделирования, рассчитаны основные функциональные 
параметры предлагаемой конструкции. 

Описание испытательного стенда. Стенд для ресурсных испытаний поршневых гидроцилиндров с 
рекуперацией энергии разрабатывался на основе ранее предложенного метода испытаний с рекуперацией энер- 
гии объемных гидравлических машин вращательного действия [3, 4]. Данное решение подразумевает, что гид- 
равлический мотор через систему механического привода возвращает энергию на вал гидравлического насоса. 
Испытываемые гидроцилиндры размещаются в напорной магистрали между гидронасосом и гидромото- 
ром [5, 6]. Такие испытания позволяют существенно уменьшить сам стенд и сэкономить значительное количе- 
ство энергии за счет ее рекуперации. 

Принципиальная гидравлическая схема стенда для испытаний поршневых гидроцилиндров изображена 
на рис. 1. 
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Рис. 1. Принципиальная гидравлическая схема стенда для ресурсных испытаний поршневых гидроцилиндров 
с рекуперацией энергии 


На рис. | показано, что в напорную гидролинию между гидронасосом и гидромотором устанавливают- 
ся испытываемые гидравлические цилиндры Ц] и Ц2, связанные между собой механической передачей МПЗ. 

Стенд функционирует следующим образом. Электрический двигатель ЭД посредством механической 
передачи МПТ приводит во вращение вал гидравлического насоса Н. Энергия, сообщаемая при этом рабочей 
жидкости, переносится ею по гидролинии |—2 на вход гидравлического распределителя Р, который направляет 
ее, например, по магистрали 4—5—6 в поршневую полость гидравлического цилиндра Ц1. 

Гидроцилиндр Ц1 посредством механической передачи МПЗ передает полученную от рабочей жидко- 
сти энергию на шток гидравлического цилиндра Ц2, который в этом случае выполняет функцию насоса и пере- 
дает энергию рабочей жидкости, находящейся в его поршневой полости. 

Из поршневой полости гидроцилиндра Ц2 рабочая жидкость по гидролинии 12—13—14—15—16 подается 
на вход гидравлического мотора М, который преобразует полученную от рабочей жидкости энергию в энергию 
вращения вала. Вращение вала гидромотора М посредством механической передачи МП2 передается на вал 
гидронасоса Н. Передача МП2 рассчитана таким образом, что частота вращения, передаваемая на вал гидрона- 
соса Н от вала гидромотора М, несколько выше той частоты, с которой первичный источник энергии (электро- 
двигатель ЭД) вращает вал гидронасоса. Это способствует затормаживанию вращения вала гидромотора М. Как 
следствие, повышается давление на входе и открывается предохранительный клапан КП2. Давление возрастает 
также в гидроцилиндрах Ц] и Ц2, что обусловливает их работу в соответствующем режиме. 

Когда шток гидроцилиндра Ц1 выдвигается на полный ход, на распределитель Р подается команда о 
смене позиции — и рабочая жидкость, поступающая на вход распределителя Р от гидронасоса Н, направляется 
по гидролинии 9-—8—7 в штоковую полость гидроцилиндра Ц1. Это вызывает обратное движение его поршня, но 
система рекуперации энергии функционирует, как и при прямом ходе поршня. 

Результаты исследования 

Математическое моделирование гидравлической системы стенда. Разработаем математическую 
модель предлагаемой рекуперативной системы испытаний поршневых гидравлических цилиндров. В качестве 
основы используем теорию объемной жесткости [11-13] с учетом приведенных коэффициентов объемной 
жесткости гидравлических элементов. Такой подход позволяет более точно моделировать систему, приближен- 
ную к реальным условиям эксплуатации гидроцилиндров [14—19]. Особое внимание при моделировании гид- 
равлических приводов следует уделять определению приведенного коэффициента объемной жесткости гидрав- 
лических магистралей. Его величина для металлических трубопроводов рассчитывается по известным зависи- 
мостям, а для рукавов высокого давления (РВД) определяется экспериментально [20, 21]. 

В соответствии с теорией объемной жесткости уравнение приращения давления в любой точке гидрав- 
лической системы может быть определено по уравнению 
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а =. () 0,м- У Ом) а, 


где Ух Оки У, Озых — суммарные входящие и выходящие из рассматриваемого (1-го) объема системы расходы 





рабочей жидкости за время 4; С; — приведенный коэффициент объемной жесткости выбранной области гидро- 
системы. 

Условно разделим гидравлическую систему стенда (рис. 1) узловыми точками. Примем как точку 0 
объем гидравлического бака с давлением, равным атмосферному, и запишем уравнения для определения давле- 
ния в выбранных узловых точках. 


ар: = С1(9н- 91-2 — О1-з) Е, 
ар» = С›(01-2 — 9—4 ) а, 
арз = Сз(@1-з — Окпа АЕ, 
арл = С4(02-4 — 94-5 ДЕ, 
ар = С5(94-5 — 95-в )4Е, 
аРь = Си1 (95-6 — "1 АЕ, 

ар = Сита (а Лшт — 97-в ) АЕ, 
@рз = Сзв(97-в — Оз-о)аЕ, 
Яро = Со(@з-9 — 99-10), 

АРло = Сто (9-10 — 910-11), 

@Р11 = С11(910-11 — 911-18 — ока) а, 
@р17 = С17(9окла — 917-20 — @ок2) 4, 
@Р2о = С20(@917-20 — 920-19), 

АРло = Сит2 (020-19 — Ип2Лишт) АЕ, 
Ар12 = Си2 (2, — 912-13 Е, 

АРлз = С1з(012-13 — 913-14), 

@р14 = С14 (13-14 + @окз — @ок1 Е, 
АР15 = С15(@ока + Фок2 — 915-16 )АЕ, 
АРлв = С16 (15-16 + @ок5 — @кп2 — @м Е, 
4рав = Сиз(011-1в + Ом — Одь — Оокз — Оокз) а. 


Здесь ар1...4р5 — приращения давления в характерных точках напорной гидролинии цилиндра Ц] за время (4; 
Ярт...Ар1о — приращения давления в характерных точках сливной гидролинии цилиндра Ц1 за время 41; 
Яи...Ар1з — приращения давления в характерных точках гидравлической системы рекуперации энергии, вклю- 
чающей гидроцилиндр Ц2 и гидромотор М, за время 4; арб и арт — приращения давления в поршневой и што- 
ковой полостях гидроцилиндра Ц] за время 4; ар1› и ар1эо — приращения давления в поршневой и штоковой 
полостях гидроцилиндра Ц? за время 41; ар›о — приращение давления на выходе из штоковой полости гидро- 
цилиндра Ц2 за время (1; Ст... С5, Сз...Сиь Сзз...С1в И С>о — приведенные коэффициенты объемной жесткости в 
характерных точках гидравлической системы стенда; Си и С‚› — приведенные коэффициенты объемной жест- 
кости поршневых полостей гидравлических цилиндров Ц] и Ц2 [11-13]; Сши и Си — приведенные коэффици- 
енты объемной жесткости штоковых полостей гидравлических цилиндров Ц] и Ц?2 [11-13]; Он — производи- 
тельность гидравлического насоса Н; Ом — расход рабочей жидкости через гидравлический мотор М; Оок‹... 
Оок5 — расходы рабочей жидкости через обратные клапаны ОКТ1...ОК5; Окт: и Окп› — расходы рабочей жид- 
кости через предохранительные клапаны КПП и КП2; О! з, О! ›, О 4, Ол 5, Оз в, Оль Оз о, Ооо, Ою п, Чи 1, 
Ото 20, Оз 14 — расходы жидкости на соответствующих участках гидравлической системы стенда; 911 И 12 — 
скорости движения поршней гидроцилиндров Ц] и Ц2 соответственно; /, — площади поршней испытываемых 
гидравлических цилиндров Ц] и Ц2; /, иг — площади поршней испытываемых гидравлических цилиндров Ц] и 
Ц2 со стороны штоковых полостей. 
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Значения расходов рабочей жидкости, необходимые для расчета приращения давлений, определяются 
по формуле: 


2 
0; = ВУ р — Р:+1| ' $181 (ф; — Рича). 





Здесь р; и р;.1 — величины давлений на входе и выходе гидравлических сопротивлений; И — площадь живого 


сечения соответствующего сопротивления; и — коэффициент расхода сопротивления; р — плотность жидко- 
сти. 

Для участков гидравлических линий (линейных гидравлических сопротивлений) приведенный коэффи- 
циент расхода определяется по формуле: 


1 
и = |еИ = ] , 
Ры = 
1 
где 4 и [ — внутренний диаметр и длина соответствующего участка трубопровода; ХЛ, — величина 


коэффициента гидравлического трения участка трубопровода, определяемая с учетом режима течения рабочей 
жидкости и свойств трубопровода. 

Приведенные коэффициенты объемной жесткости металлических трубопроводов определяются по 
формуле [11-13]: 








__4_ _ ЕЯ 
1 — 2 ЕЕ > 
па71 а“ 7 
1+5. 


где 4 и / — диаметр рассматриваемого трубопровода и его длина; д — толщина стенки трубы; Ел и Е; — вели- 
чины модуля упругости жидкости и материала стенки трубопровода. 

Величину приведенного коэффициента объемной жесткости РВД и трубопроводов, изготовленных из 
упругих материалов, необходимо определять экспериментально [20, 21]. 


Производительность насоса определяется с учетом его объемного коэффициента полезного действия 


Чнбн 
Он = 21 По, 


где ан — рабочий объем гидронасоса; он — частота вращения вала гидронасоса; о — величина мгновенного 





значения объемного КПД насоса. 


Рн 
По =1- (1 бояр) о 


пот 
Здесь Поло, — Номинальная величина объемного КПД насоса (принимается равной объемному КПД при номи- 


нальном давлении насоса); рт — номинальное рабочее давление гидронасоса; ри — текущая величина давле- 
ния на выходе насоса (давление в точке 1 гидравлической системы). 

Моделирование движения подвижных элементов гидравлической системы стенда. Величины ра- 
бочих зазоров обратных клапанов определяются из уравнения движения их затворов: 


Чл 1 [а ЯЙкл 


= — — — Е й — = т 
АЕ Тк 4. (Ракл Ракл) пр АЕ Ик 


где ух —_ скорость перемещения затвора клапана, тк, —— приведенная масса затвора клапана, Ик —_ Перемеще- 








ние (величина рабочего зазора) затвора клапана; 4„„ — диаметр отверстия клапана; Ён, — сила воздействия 
пружины на затвор клапана; { — время. 
Закон движения роторов гидронасоса Н и гидромотора М описывается уравнением их движения: 


м 1 
а = Г [мм (ам — Р2м) — Мм], 
он 1 


—— = — (Мэдып1 + Ммып2 — Ун(Ф1 — Ра). 
ЧЕ ]р 


Здесь им, ин — характерные объемы мотора М и насоса Н соответственно; хм и ®хн — угловые скорости вра- 
щения валов гидромотора М и насоса Н; Л, Ли — центральные моменты инерции роторов гидронасоса и гидро- 
мотора; р1м — давление на входе гидравлического мотора М; р›м — давление на выходе гидромотора М; р. — 
давление на входе насоса Н (давление в баке, принимается равным атмосферному давлению); р! — давление на 
выходе гидравлического насоса Н; Мэд и Мм — крутящие моменты, создаваемые электродвигателем и гидро- 
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мотором соответственно; 1мп! И &мп2 — передаточные отношения механических передач МП? (от вала гидромо- 
тора на вал гидронасоса Н) и МП] (от вала электрического двигателя на вал гидронасоса Н) соответственно. 

Механическая передача МП2 обеспечивает соотношение скоростей вращения валов гидромотора М и 
гидронасоса Н, описываемое формулой: 

«м = @н@мп2. 

Моделирование механической схемы стенда. Для определения величин скоростей перемещения 
плунжеров гидравлических цилиндров рассмотрим кинематическую схему механической передачи стенда 
МПЗ (рис. 2). 


{ 








Рис. 2. Кинематическая схема механизма передачи движения между испытуемыми цилиндрами 


Передача работает следующим образом. Корпуса гидроцилиндров Ц1 (Ё1) и Ц2 ([.2) шарнирно крепят- 
ся соответственно в точках Си (2, а их плунжеры шарнирно соединяются в точках В] и В2 с коромыслом, 
имеющим ось вращения в точке А. 

Допустим, что ведущий гидроцилиндр (гидравлический двигатель) является гидроцилиндром [Гщт, а ра- 
бочий цилиндр (гидравлический насос) является гидроцилиндром Г. Скорость плунжера гидроцилиндра Ги 
устанавливается соответственно расходу Оз 5 рабочей жидкости, поступающей в его поршневую полость. Тогда 
скорость движения его плунжера можно определить по формуле 


к 956 
1 — . 
° Л 


Движение плунжера гидроцилиндра Гл! через шарнир В! передается на коромысло АД, которое враща- 





ется вокруг точки А. Разложив скорость соединения В! на радиальную 17. и тангенциальную 7.1, определим 
значение тангенциальной составляющей: 
7-1 — Ил " шву Е 


Тогда угловую скорость ©_р вращения коромысла А) можно определить из выражения 
91 


®др = 
в 


где 7, — длина радиуса, соединяющего точку А вращения коромысла с шарниром ВБ\. 

Через коромысло АР движение плунжера гидроцилиндра Г, передается плунжеру гидроцилиндра Г», 
связанного с коромыслом АД посредством шарнира Б>. В этом случае тангенциальную скорость шарнира В» 
определим по формуле: 

Ит2 — Т2 ФАР, 
где 7› — длина радиуса, соединяющего точку А вращения коромысла с шарниром ВБ. 

Спроецировав величину тангенциальной скорости движения шарнира В. на направление перемещения 

плунжера гидроцилиндра /..2, определим скорость его перемещения 





пе”. 
5тВ> 
Величины углов В: и В› определятся по закону косинусов из треугольников АСВ! и АС>В» соответ- 
ственно: 
Сов, — Е Е, | 
2 АТ 
а [ща + т? — [2 
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Здесь [11 И [2 — расстояния между осями крепления соответствующих гидроцилиндров к раме стенда 
и к коромыслу с учетом степени выдвижения их плунжеров; [1 и [> — расстояния между осью вращения коро- 
мысла и точками крепления к раме стенда шарниров корпусов гидравлических цилиндров. 

Обсуждение и заключения. Представленная в статье система уравнений показывает, каким образом 
определяется приращение давления в выбранных узловых точках системы рекуперации энергии (в частности, 
как приращение зависит от времени, приведенного коэффициента объемной жесткости, расхода рабочей жид- 
кости, площади поршней). Величины скоростей перемещения плунжеров гидравлических цилиндров определе- 
ны согласно кинематической схеме механической передачи стенда. Итоги исследования позволяют утверждать, 
что, благодаря представленному в статье решению, результаты ресурсных испытаний гидравлических цилин- 
дров будут адекватно отражать их работу при номинальном режиме эксплуатации. 
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МАШИНОСТРОЕНИЕ И МАШИНОВЕДЕНИЕ 
МАСШИЛЕ ВОПОШУС АХО МАСНГУЕ ЭСТЕМСЕ 





УДК 620.179 №рз://401.0г©/10.23947/1992-5980-2019-19-3-250-255 





Обоснование выбора способов контроля утечек на газопроводе" 


А. Д. Лебедь’, С. П. Глушко?” 


'? Кубанский государственный технологический университет, г. Краснодар, Российская Федерация 
З@аесйоп гайопае Гог 1еаКасе топйогте ш са$ рреште””" 


А. О. Т.еьеа', $. Р. СТа$В Ко?” 


7 Кибап Э{ае Тесвпооз1са! Ошуегзйу, Ктазподаг, Вазз1ап Еедеганоп 


Введение. Анализируются эффективные методы обнару- птодиспоп. ЕЁсет{ [еак деесноп те@о4$ ап газ Ном 
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Введение. В связи с высокими темпами развития газовой промышленности, ростом удельного веса при- 
родного газа относительно других источников энергии на внутреннем российском рынке и увеличением потреб- 
ления газа на мировом рынке актуальными задачами являются снижение потерь газа и повышение безопасности 
эксплуатации газопроводов [1]. Поэтому важнейшая проблема заключается в сохранении исправного состояния 
линейных участков промысловых и магистральных газопроводов (МГ). Подземные газопроводы, эксплуатируе- 
мые при нормальных режимах, могут сохраняться несколько десятилетий в рабочем состоянии. Сохранению ис- 
правного состояния подземных и надземных МГ способствуют контроль их технического состояния и своевре- 
менное проведение регламентных работ. 

Дефекты на МГ появляются в результате коррозии и, реже, из-за механических повреждений. Определе- 
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ние мест коррозии и повреждений связано с временными и материальными затратами. Вскрытие газопровода для 
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его непосредственного визуального обследования экономически невыгодно или невозможно. Обследование лишь 
наружной поверхности трубы обычно не дает желаемых результатов. Поэтому актуальным является мониторинг 
состояния подземных и надземных промысловых, МГ без их вскрытия. Решение этой задачи осложняется значи- 
тельными техническими трудностями, однако современные методы и средства измерительной техники позволяют 
их преодолеть. Эти средства различаются по следующим признакам и параметрам: 

® по физическим явлениям, заложенным в основу работы; 

е по принципу действия; 

е по чувствительности; 

® по области применения; 

® по локальности или глобальности диагностирования. 

В системах обнаружения утечек газа в трубопроводах используются органолептические и инструмен- 
тальные методы. Самый простой метод — это фиксация запаха. При визуальных методах контроля применяют 
вспениваемую мыльную эмульсию, фиксируют появление наледи или снежной шубы, появление желтой травы 
летом или бурого снега зимой. Безусловно, инструментальные методы являются более совершенными и точными. 
В их основе лежит применение приборов-газоанализаторов автоматического и ручного действия. В качестве ав- 
томатических газоанализаторов широко применяются датчики загазованности — это высокоточные измеритель- 
ные устройства, используемые при непрерывном автоматическом контроле количества газов, содержащихся в 
воздухе. Сигнализаторы загазованности применяются в промышленных, бытовых и коммунальных помещениях. 
Датчики газа выпускаются в виде отдельных приборов, которые могут размещаться автономно и реагировать на 
изменение концентрации определенных газов. Они мгновенно реагируют на повышение концентрации газов, ис- 
ходящих от газотранспортной системы. 

Современные модели датчиков оснащены жидкокристаллическими дисплеями, на которых отображаются 
измеряемые показатели. Встроенная память позволяет фиксировать полученные данные и хранить все результаты 
измерений. Также на дисплей выводится информация о работоспособности прибора и имеющихся неполадках. 
Эта функция позволяет своевременно менять изношенные элементы. 

Датчики контроля загазованности делятся на два основных вида: бытовые и промышленные. Сигнализа- 
торы бытового применения обеспечивают ряд реакций на превышение установленной концентрации газа. Основ- 
ные задачи промышленного датчика — измерение и отображение на дисплее показателей загазованности. Про- 
мышленные приборы применяются для решения более значимых задач и используются в составе автоматизиро- 
ванных систем вместе с пультом управления, на котором отображена информация, поступающая от датчиков газа. 

По виду контролируемого газа различают метановые, пропановые сигнализаторы, а также сигнализаторы, 
предназначенные для измерения концентрации угарного газа или окиси углерода. Наряду с этим существуют ком- 
бинированные (многокомпонентные) сигнализаторы. Газовые датчики могут быть стационарными, которые рабо- 
тают только от электросети и располагаются вблизи от источника электропитания, или переносными, питание 
которых осуществляется от встроенных аккумуляторов. 

Датчик загазованности является многофункциональным устройством и имеет функции звукового и свето- 
вого оповещения. Например, в системе САКЗ-МК-1, если прибор отреагировал на загазованность, у него загорает- 
ся индикатор и срабатывает громкий звуковой сигнал (рис. 1). Датчики загазованности имеют релейный выход 
для подключения дополнительных электрических устройств: оповещателей, вентиляторов, механизмов и прибо- 
ров пожарного и диспетчерского пультов. 







На 
дублирующие 
сигнализаторы 


“220 В 


КЗГЭМ-У Кран 


Направление 
подачи 
среды 


=> 


Рис. |. Система сигнализации загазованности САКЗ-МК-1 


Машиностроение и машиноведение 


= 


БИр://уезбиК.Чоп$ва.га 


р 


Лебедь А. Д. и др. Обоснование выбора способов контроля утечек на газопроводе 


Гере4 А. О., её. 41. з@еспоп гапопае ог [еаКаге топйотп? т газ ррейпе 





Датчики загазованности, согласно СТО Газпром 5.37-2011, используются в узлах измерений расхода и 
количества горючего природного газа, поставляемого от объектов промысловой подготовки, подземных храни- 
лищ, газоперерабатывающих заводов в магистральные газопроводы и далее транспортируемого потребителям 
Российской Федерации, в страны ближнего зарубежья, утилизируемого и потребляемого на собственные техно- 
логические нужды. На всем протяжении магистральных трубопроводов датчики загазованности использовать 
невыгодно. Кроме того, в труднодоступных местах расположения магистральных трубопроводов применение 
органолептических и инструментальных методов обнаружения утечек газа затруднительно. 

Материалы и методы. Помимо вышеперечисленных методов обнаружения утечек газа в трубопрово- 
дах используют следующие [2]: 

® метод анализа профиля давления; 

е объемно-балансовый метод; 

® метод анализа акустической эмиссии; 

® метод переменного перепада давления; 

® механический метод с применением тахометрических расходомеров; 

® ультразвуковой метод измерения расхода. 

Согласно методу выявления утечек по анализу профиля давления, проводятся моделирование распре- 
деления давления вдоль трубопровода и статистический анализ этого распределения. При утечке газа расход на 
контролируемом участке становится больше первоначального расхода на этом участке, а расход на участке по- 
сле места появления утечки становится меньше первоначального расхода. Вследствие этого перепад давления 
на участке до места утечки увеличивается, а после утечки уменьшается, что приводит к появлению излома в 
приведенном профиле давления. Метод выявления утечек по анализу профиля давления работает только в ста- 
ционарном режиме, так как многие факторы вызывают похожие изменения. Существенным недостатком данно- 
го метода является низкая точность и наличие ложных срабатываний. 

Для уменьшения ложных срабатываний отклонения должны быть зафиксированы как минимум в двух 
соседних точках. Для этих же целей используется усредненный профиль распределения давления, который яв- 
ляется квазистационарным профилем. Этот профиль получается в результате фильтрации давления в точках 
трубопровода. Усредненный профиль давления не является постоянным. Он изменяется, но медленнее, чем ре- 
альное давление, что приводит к увеличению времени обнаружения утечек. Данный метод из-за неспособности 
обнаружения небольших утечек и большой погрешности при определении координаты места утечки в настоя- 
щее время в основном используется только совместно с другими методами. 

Объемно-балансовый метод контроля основан на том, что при образовании утечки расход на входе ста- 
новится больше расхода на выходе. Кроме того, этот метод учитывает количество газа в самой трубе, которое 
при появлении утечки уменьшается. Для реализации этого метода необходимо измерять расход на границах 
контролируемого участка с помощью высокоточных приборов. 

Параметром контроля при объемно-балансовом методе является не давление газа, а его нормализован- 
ный расход. Контроль участка газопровода осуществляется путем определения разности нормализованных обЪ- 
емов газа, входящих и выходящих из участка между двумя локальными расходомерами. Этот метод позволяет 
диагностировать как быстро развивающиеся разрывы в трубе, так и медленно развивающиеся утечки, а также 
утечки на очень больших участках трубопровода между расходомерами. Минимальная величина диагностиру- 
емой утечки определяется, в первую очередь, погрешностью измерения расхода и, при имеющихся в настоящее 
время средствах измерения расхода, находится на уровне 0,5—1,0 %. Точность контроля утечек зависит от точ- 
ности расходомеров. Невозможно обнаружить утечку, расход которой меньше, чем погрешность измерения. В 
данном случае имеет значения не абсолютная погрешность, а только погрешность одного прибора относитель- 
но другого. 

Объемно-балансовый метод не позволяет определять координату утечки — это его существенный не- 
достаток. Кроме того, расход объемного баланса по участкам позволяет определять возможное место утечки в 
газопроводе только с точностью оценки этого параметра до конкретного участка. Это, в свою очередь, приво- 
дит к дополнительным затратам времени и потерям газа с момента утечки до момента её устранения. 

Метод анализа акустической эмиссии лидирует в индустрии автоматического обнаружения утечек газа 
в трубопроводах [3]. Этот метод основан на регистрации и анализе акустических волн, возникающих при пла- 
стическом деформировании и микроразрушении материала трубопровода, а также при истечении газа через 
сквозные отверстия в трубе. Для приема сигналов акустической эмиссии применяются пьезоэлектрические пре- 
образователи и быстродействующие измерители давления. 
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Преимущества данного метода: 

® высокая чувствительность к развивающимся дефектам; 

® малое время обнаружения; 

® высокая точность определения координат места утечки. 

Недостатком метода является трудность выделения сигнала акустической эмиссии на фоне шумов и 
помех. Для повышения помехоустойчивости и уменьшения количества ложных срабатываний применяются 
специальные методы обработки принимаемых сигналов. 

При разработке системы обнаружения утечек газа на магистральных газопроводах и узлах учета расхо- 
да газа на технологическом объекте должны решаться следующие задачи: 

е достижение независимости результатов измерения от изменения параметров среды за счет контроля 
давления газа; 

е повышение точности контроля утечек и расхода газа. 

Измерение расхода газа является важной задачей в газовой промышленности, так как система учета ко- 
личества газа невозможна без средств измерения расхода. На объектах газовой промышленности для измерения 
расхода газа и его количества применяются тахометрические расходомеры и счетчики. В последнее время в 
индустрию контроля газопроводов также активно внедряются наиболее прогрессивные ультразвуковые расхо- 
домеры [4]. Однако в настоящее время для измерения расхода газа наибольшее распространение получил метод 
переменного перепада давления на сужающем устройстве (СУ). Суть метода состоит в измерении перепада 
давления до и после СУ, установленного в потоке газа. Объемный расход газа О для специальных СУ опреде- 
ляют по формуле: 


О = СЕМ е- : 
р 
где С — коэффициент истечения, представляющий собой отношение действительного расхода к теоретическо- 
му; Е — коэффициент скорости входа; = — коэффициент расширения, учитывающий увеличение удельного 
объема газа; Ру — площадь СУ, м. Ар — перепад, создаваемый СУ, Па; р — плотность газа, кг/м". 

Достоинства метода переменного перепада давления: 

1. Универсальность. Данный метод применяется для измерения расхода практически любых сред: 
жидкости, газа, пара. Для вязких жидкостей применяются СУ специальной формы. 

2. Низкая первоначальная стоимость. Стоимость расходомера, основанного на методе переменного 
перепада давления, складывается, исходя из стоимости СУ, импульсных линий и датчика дифференциального 
давления. 

3. Простая методика поверки. Для периодической поверки расходомеров требуется измерить геомет- 
рические размеры СУ и поверить датчик дифференциального давления. 

4. Отсутствие движущихся частей. 

5. Возможность измерения расхода в условиях высокого давления. Давление в трубопроводе может 
достигать 40 МПа. 

6. Возможность измерения расхода в широком диапазоне температур измеряемой среды — от минус 
200°С до плюс 1000°С. 

7. Широкий диапазон типоразмеров. Рассматриваемый метод используется на трубопроводах в диапа- 
зоне диаметров условного прохода от 15 до 2000 мм. 

Недостатки измерения расхода газа методом переменного перепада давления на СУ: 

1. Узкий динамический диапазон. Стандартный динамический диапазон СУ приблизительно 1:3. Та- 
кое ограничение связано, в первую очередь, с квадратичной зависимостью между расходом и перепадом давле- 
ния на СУ. Увеличить динамический диапазон позволяет использование высокоточных датчиков дифференци- 
ального давления. 

2. Высокая стоимость эксплуатации. Расходомеры, основанные на СУ, требуют периодического об- 
служивания, которое включает в себя: измерение геометрических размеров сужающего устройства, прочистку 
импульсных линий, прогрев импульсных линий, установку нуля на датчике дифференциального давления. 

3. Низкая точность измерений. Погрешность измерений обычно составляет 3,0-3,5 %. 

Принцип действия тахометрических расходомеров основан на зависимости скорости вращения преоб- 
разователя, установленного в трубопроводе, от расхода газа. У расходомеров типа «счетчик» вал преобразова- 
теля расхода через редуктор связан со счетным механизмом, что позволяет измерять количество прошедшего 
газа. Достоинствами таких приборов являются быстродействие и большой диапазон измерений. 


Машино строение и машиноведение 


й 5 


БИр://уезбиК.Доп$ва.га 


254 


Лебедь А. Д. и др. Обоснование выбора способов контроля утечек на газопроводе 


Гере4 А. О., её. 41. З@еспоп гапопае [ог ([еаКаге топйотт?г т гах ррейпте 





Погрешность счетчиков составляет 0,5-1,5 %, использование тахометрических преобразователей поз- 
воляет снизить погрешность преобразования расхода в частоту вращения преобразователя до 0,3 %. Основные 
недостатки тахометрических расходомеров — износ опор и наличие подвижных элементов. Недостаток по от- 
ношению к расходомерам с СУ — необходимость в поверочных установках. 

Метод, на основе которого работают ультразвуковые расходомеры, основан на явлении смещения зву- 
кового колебания движущейся средой. Этот метод имеет следующие недостатки: 

е зависимость собственной скорости ультразвуковых колебаний от физико-химических свойств изме- 
ряемой среды; 

® скорость потока усредняется вдоль ультразвукового пучка, а не по сечению трубы. 

Усреднение скорости потока заставляет разработчиков снабжать конструкцию дополнительными дат- 
чиками или отражателями, что делает расходомер более сложным. При этом возрастает вероятность ошибки 
при выходе из строя датчиков системы. 

Достоинства данного метода: 

е отсутствие перепада давления; 

® высокое быстродействие; 

® отсутствие подвижных элементов. 

Результаты исследования. Анализ показывает, что задача обнаружения и локализации (точного опре- 
деления места) утечки газа не может быть решена с использованием одного из рассмотренных методов. С уче- 
том достоинств и недостатков методов, рассмотренных в данной статье, предлагается решить задачу с приме- 
нением трех методов и поэтапно: 

® определить участок с утечкой газа между станциями методом расхода объемного баланса; 

® методом анализа акустической эмиссии установить локальный участок повреждения трубопровода 
и утечки газа; 

® с помощью ультразвуковой дефектоскопии определить координаты утечки в пределах установлен- 
ного локального участка повреждения. 

Для контроля горизонтальных и прямых участков газопровода и точной фиксации мест утечек с помо- 
щью ультразвукового сканера возможно применение различных роботов [5—7], например, видеокроулеров, ко- 
торые представляют собой мобильные и универсальные средства инспекции трубопроводов. Они подходят для 
различного применения благодаря модульной конструкции и способности осматривать внутренности труб с 
диаметрами от 150 до 900 мм. На рис. 2 показан один из таких роботов — видеокроулер Коууег 600. Дистанци- 
онное управление робота дает возможность контролировать фокус и освещение, а также направлять ход в обла- 
сти препятствий или ответвлений. 





Рис. 2. Видеокроулер Коууег 600 


Обсуждение и заключения. Авторы предлагают для повышения надежности обнаружения и локализа- 
ции утечки газа применять комбинированную трехэтапную методику, включающую последовательное исполь- 
зование метода расхода объемного баланса, метода анализа акустической эмиссии и ультразвуковую дефекто- 
скопию. В определенных ситуациях рекомендуется применение специальных роботов. Авторы считают, что 
роботы должны снабжаться датчиком обнаружения 90-градусного изгиба, установленным на головном звене, а 
также датчиками скорости и пройденного пути для обеспечения точного позиционирования. Это позволит 
уменьшить время диагностирования газопроводов различных типов и сократить затраты на контроль их состо- 
яНИЯ. 
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Анализ информации по авариям в газовом комплексе имеет большое значение, прежде всего, для оцен- 
ки показателей риска и надежности работы газотранспортной системы. Результаты анализа требуются для 
обоснования решений при проектировании объектов системы газоснабжения, а также при планировании меро- 
приятий по повышению их надежности и безопасности в процессе эксплуатации. 
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Введение. Обмолот представляет собой многофакторный 
процесс, на который влияют, в частности, скорость подачи 
растительной массы, удельный вес обмолоченного зерна и 
недомолота, показатели сепарации, засоренность и влаж- 
ность растительной массы. С этой точки зрения актуальны 
и вопросы оптимизации профиля среза транспортирующе- 
го устройства, которая позволяет проводить обмолот и 
сепарацию с наименьшими усилиями. Опытным путем 
установлено также, что на процесс обмолота влияют объем 
подаваемого зернового материала за секунду (секундная 
подача) и скорость движения агрегата. 

Материалы и методы. Исследования проводились на ис- 
пытательном стенде, оснащенном молотильно- 
сепарирующим устройством в форме однополостного ги- 
перболоида. При выполнении представленной работы ши- 
рина барабана была разделена на три зоны, а длина — на 
пять ячеек. Зерно пшеницы, полученное при обмолоте и 
сепарации в каждой зоне и ячейке, поступало в отдельные 
контейнеры. Солома собиралась отдельно. Затем зерно и 
солома взвешивались. Полученные в итоге данные обраба- 
тывались статистическими и математическими методами. 
Результаты исследования. Рассчитаны длины каждой 
ячейки исходя из отношения суммарной сепарации и коли- 
чества зерновой массы, приходящейся на единицу длины. 
Вычислены длины каждой ячейки в зависимости от изме- 
нения влажности. Показатели соотношения влажности 
массы и длин ячеек представлены в виде таблицы. По за- 
данным табличным значениям для различных уровней 
влажности зерновой массы построены графики, каждый из 
которых описан математической моделью, учитывающей 
длину и ширину барабана. Представлен усредненный про- 
филь среза для обрабатываемой растительной массы влаж- 
ностью 8%, 12%, 16%, 20%, 24%. 

Обсуждение и заключения. Анализ данных этой и более 
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ранних работ позволил сравнить графическое выражение 
зависимостей профиля среза транспортирующего устрой- 
ства от секундного объема и влажности поступающего на 
него зернового вороха. Установлено, что кривые профиля 
среза идентичны по всей длине барабана. Коэффициент 
достоверности близок к 1, что говорит о точности модели. 
Очевидна идентичность усредненных профилей среза в 
зависимости от влажности растительной массы и от се- 
кундной подачи. 


Ключевые слова: молотильно-сепарирующее устройство, 
зерновой ворох, колосовая фракция, подача растительной 
массы, транспортирующее устройство, профиль среза 
транспортирующего устройства, недомолот, сепарация, 
засоренность, влажность. 
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Введение. Улучшение качественных показателей зерна — один из важных факторов повышения уро- 
жайности [1, 2]. В процессе зерноуборочных работ, особенно обмолота, происходит сильное повреждение зе- 
рен [3]. В зерноуборочном комбайне зерно отделяется от соломенных примесей. На этом этапе зерна получают 
повреждения, ухудшается их качество, происходят количественные потери [4]. 

Можно сократить негативное воздействие названных факторов, если в процессе обмолота увеличить 
путь прохождения массы равномерным слоем по всей длине барабана. Проведенные ранее исследования пока- 
зали, что обмолачиваемое зерно проходит более плотным слоем в середине барабана [5], и в этой зоне увеличи- 
вается сепарация. Ближе к краям барабана масса значительно менее плотная, и объем отделенного зерна здесь 
в 2—3 раза меньше, чем в середине [6, 7]. С этой точки зрения наиболее перспективной представляется кон- 
струкция молотильно-сепарирующего устройства тангенциально-аксиального типа, который равномерно рас- 
пределяет обмолачиваемую массу по всей длине барабана. 

Потери при обмолоте и сепарации зависят от качественных показателей зерна и от технических харак- 
теристик барабана [8, 9]. Актуальна также оптимизация профиля среза транспортирующего устройства, которая 
позволяет проводить обмолот и сепарацию с наименьшими усилиями. Профиль среза транспортирующего 
устройства — это часть транспортной доски сепаратора, предназначенная для равномерного и непрерывного 
перемещения колосовой фракции по ширине тангенциально-аксиального молотильно-сепарирующего устрой- 
ства в виде гиперболоида вращения. 

В свое время была установлена зависимость профиля среза от частоты вращения режущего рабочего 
органа, высоты его установки и угла наклона режущих элементов жатки [10]. Кроме того, опытным путем уста- 
новлено, что на процесс обмолота также влияют объем подаваемого зернового материала за секунду (секундная 
подача) [8], а также скорость движения агрегата [11]. 

Цель исследования — создание математической модели процесса обмолота при равномерной подаче 
обрабатываемой массы по всей длине молотильно-сепарирующего барабана. 

Материалы и методы. Исследования проводились на испытательном стенде, оснащенном молотиль- 


но-сепарирующим устройством в форме однополостного гиперболоида (рис. 1). 
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Рис. 1. Схема стенда, оснащенного молотильно-сепарирующим устройством: вид сбоку (а); вид сверху (6) 


Стенд состоит из электродвигателя 1. От него через ременную передачу на приводной барабан 4 пере- 
дается вращение. На барабане установлено молотильно-сепарирующее устройство. Его форма представляет 
собой вращающийся однополостный гиперболоид 7. Молотильно-сепарирующее устройство 7 состоит из от- 
бойных битеров 6 и деки 5. Испытательный стенд защищен кожухом 8. В загрузочную зону 2 поступает коло- 
совая фракция. Эта масса по транспортерной ленте 3 подается к молотильно-сепарирующему устройству 7. 

При выполнении научной работы ширина барабана была разделена на три зоны, а длина — на пять яче- 
ек. Зерно пшеницы, полученное при обмолоте и сепарации в каждой зоне и ячейке, поступало в отдельные кон- 
тейнеры для сбора. Солома собиралась отдельно. Затем зерно и солома взвешивались. 

Данные, полученные в результате проведенного исследования, обрабатывались статистическими и ма- 
тематическими методами [10, 11]. 

Предварительные исследования [6, 7, 12] позволяют сделать следующее утверждение: по мере продви- 
жения обмолачиваемого зернового вороха в молотильно-сепарирующем устройстве суммарная сепарация зерна 
увеличивается от первой зоны к третьей. 

В рамках представленной работы профиль среза определен экспериментально. 

Задача исследования — сравнить зависимости профиля от: 

— объема подаваемого в секунду растительного материала; 
— влажности растительного материала. 

Ранее установлена взаимосвязь профиля среза и величины секундной подачи [13, 14], поэтому необхо- 
димо установить аналогичную зависимость профиля и уровня влажности обмолачиваемой массы. В табл. 1 
представлены абсолютные величины сепарации по ячейкам в зависимости от влажности [12]. 

Таблица 1 
Абсолютные суммарные значения сепарации по ячейкам в зависимости от влажности массы 





Результаты и обсуждение. Для определения профиля среза должно быть выполнено следующее усло- 
вие: сохранение постоянного (с01п5@) количества обрабатываемой массы на единицу длины молотильно- 
сепарирующего устройства. Величина постоянной определяется отношением суммарной сепарации к длине 
барабана Г, = 1200 мм. Методика определения была применена ранее при построении профилей среза, соответ- 
ствующих различным величинам секундной подачи [14]. 

Для расчета длины каждой ячейки необходимо определить отношение суммарной сепарации и количе- 
ства зерновой массы, приходящейся на единицу длины: 
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где ) 5 — величина абсолютной суммарной сепарации в 1-Й ячейке; соп$ — объем обрабатываемой массы, 
1 


приходящейся на единицу длины МСУ при заданной влажности растительной массы. 
Подставляя значения из табл. | в формулу (1) и вычисляя длину каждой ячейки в зависимости от изме- 
нения влажности, получаем результаты, приведенные в табл. 2. 


Таблица 2 


Соотношение влажности массы и длин ячеек 


Г 





Суммарная величина длин всех ячеек >.( равна ширине барабана В = 1000 мм. 


Значения величин длин (; откладываются на графике по абсциссе и соответствуют 1-м ячейкам длины 


барабана Г. 

По заданным табличным значениям для каждой из величин влажности зерновой массы были получены 
приближенные функции Г(х). Каждый график, представленный на рис. 2, описан математической моделью за- 
висимости длины барабана (Г.) от его ширины (В). 
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Рис. 2. Зависимость профилей среза от показателя влажности обрабатываемой массы: 
8% (а): 12% (6); 16% (в); 20% (г); 24% (9) 
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Из рис. 2 видна идентичность графиков. Следовательно, можно построить усредненный профиль среза 
для всего диапазона влажностей обрабатываемого зернового вороха. С этой целью величины длин каждой из 
ячеек суммируются и делятся на общее количество величин влажностей. По аналогии с предыдущими графи- 
ками строится усредненный профиль среза (рис. 3). 
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Рис. 3. Усредненный профиль среза для обрабатываемой растительной массы влажностью 8%, 12%, 16%, 20% и 24% 


Данные рис. 3 сравниваются с полученным ранее усредненным профилем среза, обусловленного се- 
кундной подачей (рис. 4). 
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Рис. 4. Усредненный профиль среза для обрабатываемой растительной массы в зависимости от секундной подачи 


Обсуждение и заключения 

1. Из графиков рис. 2 видно, что при изменении влажности обрабатываемой растительной массы кри- 
вые профиля среза идентичны по всей длине барабана. 

2. Коэффициент достоверности аппроксимации В показывает степень соответствия трендовой модели 
исходным данным. Коэффициент достоверности близок к 1, что говорит о точности модели. 

3. Очевидна идентичность усредненных профилей среза в зависимости от влажности растительной 
массы (рис. 3) и от секундной подачи (рис. 4). 
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Введение. Технологии математического и логического мо- 
делирования решения задач по существующей практике 
их распространения распределяются на два направления: 
широко распространенное математическое моделирование 
и информационное логическое моделирование, которое в 
настоящий момент развито недостаточно, в особенности 
для сложноорганизованных систем. Принципиальные раз- 
личия этих технологий, в частности для подготовки про- 
изводства обработкой резанием, в том, что логическое мо- 
делирование информационно и логически связано с си- 
стемами организации, а математическое — с процессами 
управления в системах организации. Логическое модели- 
рование используется для оперирования геометрическими 
объектами в технологических схемах их взаимодействия 
методами базирования, геометрического формообразова- 
ния в условиях статической, т. е. идеальной настройки со- 
ответствующих схем. Математическое моделирование ис- 
пользуется для оперирования материальными объектами в 
процессах управления их преобразованиями методами об- 
работки резанием, т.е. неидеально с учетом функцио- 
нально различных погрешностей. Между рассматривае- 
мыми системами организации и процессами управления в 
них существуют информационные и логические связи их 
органического единства, отрицающие их раздельное рас- 
смотрение. Для информационной детерминированной 
технологии решения задач высокого уровня автоматиза- 
ции разграничение понятий «математическое» и «логиче- 
ское» моделирование актуально, обладает научной новиз- 
ной и практической значимостью. 

Материалы и методы. Для характеристики свойств поня- 
тий «математическое моделирование», «логическое моде- 
лирование» и функций знаний, следующих из формулиро- 
вания этих понятий, используются принципиально раз- 
личные методы и соответствующие инструментальные 
средства. В основу разграничения рассматриваемых поня- 
тий положено разграничение технологий (методы, соот- 
ветствующие средства, алгоритмы, операции) решения 
прикладных задач какой-либо предметной области знаний. 
Результаты исследования. Понятия «логическое модели- 
рование» и «математическое моделирование» являются 
концептуальными общетеоретическими понятиями, обла- 
дающими инвариантными свойствами, которые необхо- 
димы для решения задач практики какой-либо предметной 
области. В соответствии с разграничениями рассматрива- 
емых понятий, технологии решения задач подразделяются 
на два типа: технология системной инженерии — в систе- 
мах организации информационных объектов и технология 


"Работа выполнена в инициативном порядке. 
"Е-плай: е.п.Ко!убепко@ттай.га 
""ТЬе гезеагсВ 15 Чопе ут те Нате оР Фе ш4ерепдете В &0. 


[итодисНоп. Тесвпо]ослез оР тафетайса| ап4 1021са| под- 
ейпе оЁ рго ет зоуше ассог4 ие о Фе ех1$Нпсх ргасйсе оЁ 
Бег 415пиНоп аге Фту14еа шо &\мо агеаз: ул4езргеа та- 
етайса| тодейпе ап шЮ]оз1са| тодейп® \сь 1$ ситепЙу 
ипдегдеуеоре4, езреслаПу юг зоршызисаже4 зу$етз. Еипда- 
шета| 41егепсез Бебуееп Фезе {1есбио]о21ез, ш рагасиаг 
ог фе тасбите ргергодисНоп, аге ай 1о21са| то4еПие 1$ 
ш®югтайопаПу ап4 |озлсаПу ге]{е4 {о огоап1тайоп зузетл$, 
ап тафетайса! тодеппо 1$ аззослайе4 \ий сопо] ргосез$- 
ез ш Ше огоаплхайоп зузепа5. Гос1са| тодепас 1$ изе4 фо 
орегже у оеотейлс об] ес{$ ш Фе 1есбпо]ос1са| зсВетез оЁ 
Фе пиегасйоп @точоВ Базше тефо4$, сеотейлс зварше ш 
а айс (14еа1) зеие оЁ Фе соггезропд® зсВетез. Ма@е- 
тайса! зниу]айоп 1$ изе4 о орегме шаепа| об]ес{$ ш Ше 
сопёго| ргосеззез оГ Фешг НапзЮптайопз$ ФгоиеВ сие 
шефод$, 1.е. ппре{есЙу, сопз14етие Паегосепеоцз$ еггогб. 
Вебхуееп Фе огоапхайоп зу$етз ип4ег зеа4у ап тапасе- 
тепё ргосеззез ш Фет, еге аге шЮюгтайоп апа 1о21са| ПпК$ 
ог Феи оггап1с ипйу, уурлсЬ депу Фет зерагже сопзл4егайоп. 
ш Фе шЮгпайоп даегиллзИс 1есбпо]осу Юг зоуше ргоб- 
1епл$ ога №12Ъ-[еуе! ащотайоп, ве 41$НпсНоп Бебмееп Фе 
сопсер оЁ “таетайса[” ап “оглсаР” тодепие 1$ гёе- 
уапе; и Ваз зслепиЯс поуебу ап4 ргасйса] 1етийсапсе. 
Маета!5 апа Мефосб. То свагаетте Фе ргорегаез оЁ Ше 
сопсер о? “тафетайса| то4ейиэ”, “оелса| тодейпэ” апа 
фе Кпо\]е4дэе псНоп$ геза[ те Кош ве Югищайноп о е- 
5е сопсер, апдатешаПу Чегет те о4$ ап арргорпае 
{0015 аге изеа. ТВе Ч1етепнайоп оЁР Ше сопсер ипдег соп- 
51Аеганоп 15 Базе оп Фе АШегепианоп оЁ {есбпо|оз1ез 
(пе@о4$, арргорпае 1ю01[5, а|>огИбип$, ореганоп$) Юг зо1у- 
ше аррПе4 ргоепл$ оГапу Кпо\е4ее доталп. 

Кебеагсй КезиП5. ТВе 14еаз о “]оелса! тодейп?” ап “та@- 
етайса! тодеПпе” аге сопсера| сепега1-Феогейса| поНоп$ 
УИ шуапапе ргорегаез гедитеа юг зо]уше ргасйса| ргоБ- 
1еп1$ оЁ апу аррИсайоп Чоташ. ш ассогдапсе у Фе 415- 
босНой Бебмееп Фезе сопсер5, Ше рго ет зо]уе {есбпо]- 
ос1ез аге Ч1ту14е4 шю мо бурез: зует епошеегис 1есбпо]о- 
су — ш Ше огсоаплхайоп ог шЮгтайоп обес зу$фепз, ап 
зузет зс1епсе — ш Фе тапасетепе ргосеззез оЁ мапзЮг- 





Колыбенко Е. Н. Разграничение понятий математического и логического моделирования 


Ко[уфепКо Е. №. Пл5хйтспоп Бемееп ше сопсерх о тайетайса[ апа озса то4дейп? 





системотехники — в процессах управления преобразова- 
ниями соответствующих материальных объектов. Эти 
направления должны существовать в информационной и 
логической связи их органического единства. 

Обсуждение и заключения. Автором разграничены поня- 
тия «логическое моделирование» и «математическое мо- 
делирование», что является важнейшим условием успеш- 
ного перехода к детерминированной информационной 
технологии высокого уровня автоматизации в решении 
задач практики какой-либо предметной области знаний, 
например, технической подготовки производства резани- 
ем. 


Ключевые слова: техническая подготовка производства, 
обработка резанием, системный анализ, информационная 
технология, моделирование решений, системная техноло- 
ГИЯ. 
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Введение. Основные положения системной инженерии, как творческой технологии комплексного ре- 
шения инженерных и организационно-управленческих задач, изложены в работе [1]. Отмечено, что чрезвычай- 
но высокая сложность и разнообразие крупномасштабных высокоавтоматизированных систем существенно за- 
трудняет использование исключительно средств математики для их определения. Концепция системной инже- 
нерии формировалась в рамках практики успешных разработок. Однако, данная работа вышла с изменением 
буквального перевода понятия «системная инженерия» на понятие «системотехника». Далее технология систе- 
мотехники сформировалась как самостоятельное направление знаний. Функции технологии системотехники — 
управление процессами преобразования материальных объектов в соответствующих системах организации. По- 
сле введения понятия «системотехника» системная инженерия в нашей стране, как концептуальная прорывная 
конкурентоспособная информационная технология создания крупных сложных информационных систем высо- 
кого уровня автоматизации, была в значительной степени утеряна. Развитие пошло по направлению проектиро- 
вания технологических процессов управления в системах организации функционально различными методами и 
средствами математического аппарата. 

Функционально различные аспекты подхода к технологии системной инженерии в познании базы зна- 
ний предметной области по семи уровням иерархии ее базовых объектов различных типов изложены в работе 
[2] на примере технологической подготовки производства (ТП) резанием по всему ее циклу. 

Функции технологии системной инженерии — это организация на инвариантной основе детерминиро- 
ванных систем высокого уровня автоматизации для информационного логического моделирования решения 
возможного множества задач в моделях базовых объектов знаний и базах данных их решений с выпуском про- 
ектной документации. Ознакомление с работами [3-13] позволяет отметить, что направление технологии си- 
стемной инженерии широко развито в части информационного и программного обеспечения систем обработки 
информации на компьютере, а также в функционально различном менеджменте. В технологической подготовке 
механообрабатывающего производства развитие пошло в направлении перевода традиционных знаний в элек- 
тронную форму для восприятия компьютером. На основе таких преобразований формы знаний организованы 
малоэффективные не имеющие выхода в САПР диалоговые системы решения задач практики низкого уровня 
автоматизации. Технология системотехники развита широко. 

Глубинные свойства понятий приведены в работе [14]: «...одним из свойств сознания человека являет- 
ся способность применять ранее накопленные знания для решения возникающих логических проблем»; 
«...понятие можно определить и как некоторую языковую конструкцию, имеющую определенный смысл, т. е. 
образное содержание». Там же: «Объект, с целью изучения которого проводятся исследования, называется ори- 
гиналом, а объект, исследуемый вместо оригинала для изучения определенных свойств, называется моделью». 
Там же: «Моделирование есть метод (или процесс) изучения свойств объектов-оригиналов посредством иссле- 
дования соответствующих свойств их моделей». 

Если свойства реального объекта описать словами, математической формулой, чертежом или рисун- 
ком, т. е. формализовать их в модели, то такую модель можно назвать информационной. Только знания инфор- 
мационной природы пригодны для обработки на компьютере. Системный анализ знаний по технологии систем- 
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ной инженерии труден для восприятия, но его результаты информативны и одновременно лаконичны. Исходя 
из необходимости встраивания формализованных понятий в структуру базовых объектов знаний различных ти- 
пов предметной области, содержание понятий должно достигать максимальной строгости. Принципиально 
важной является авторская актуализированная трактовка существующих и принципиально новых понятий. 

Материалы и методы. Введем актуализированное понятие системных принципов. Системные прин- 
ципы — это информационные логические утверждения подхода к решению задач исследования технологиче- 
ских систем, процессов функционально различного назначения в условиях соответственно статической, дина- 
мической настройки, проверенные множественной практикой их использования. 

Далее сформулируем системные принципы. 

1. Необходимым условием упорядоченного технологического процесса функционально различной при- 
роды (технической, физической, химической, биологической) преобразования какого-либо объекта является 
наличие системы управления. Функция управления преобразованием всегда конкретна, реализуется в отноше- 
нии соответствующей технологической системы. Из этого следует, что необходимым условием существования 
системы управления преобразованиями какого-либо объекта является наличие соответствующей технологиче- 
ской системы его организации. Сошлемся на известное и справедливое высказывание: «Нельзя хорошо управ- 
лять плохо организованным объектом». При этом понятие «организация» далее рассматриваем на основе соот- 
ветствия понятию «структура» в различных аспектах его толкования на примере технологической подготовки 
механообрабатывающего производства. Решение задач ТПП автором организовано на возможном множестве 
деталей, объектов производства, станков, режущих инструментов, приспособлений для установки режущих ин- 
струментов и объектов производства; при этом используем обозначение {...} — множество. 

Введем актуализированные понятия. 

» Структура — это строение (организация) какого-либо объекта, определенное методами декомпози- 
ции в составе элементов структуры и методами синтеза наложением связей между элементами структуры раз- 
личных уровней и одного уровня структуры, в основу которого положены параметры проектного качества ос- 
новных элементов интеграции, дезинтеграции структуры. 

е Основные элементы интеграции, дезинтеграции структуры — это элементы знаний, характери- 
зующие функционально различные элементы структуры объектов ({детали}) стадии конструкторской подго- 
товки производства (КПП) и объектов ({исходные заготовки!) стадии технологической подготовки производ- 
ства (ТПО) а также связи между элементами структуры. 

е Проектное качество основных элементов интеграции, дезинтеграции в структуре объектов 
({детали}) стадии КПП и в структуре объектов ({исходные заготовки }) стадии ТПП — это функция, под- 
лежащая безусловному исполнению для совокупности параметров проектного качества, характеризующая воз- 
можное множество различных свойств объектов в их распределении по различным элементам структуры объ- 
ектов и связям между элементами структуры объектов. 

е Базовый объект знаний— это материал, который определен на основе понятия «структура» в ор- 
ганическом единстве его информационных и логических связей с понятиями «содержание» и «форма», отрица- 
ющем их раздельное рассмотрение. 

е Информационный «ведомый» (основной) объект преобразования — это первичный материал, 
информационное отображение знаний которого возможно только в единой среде одного уровня структуры тех- 
нологической схемы на основе взаимодействия элементов структуры «ведомого» объекта в общем случае с 
группой элементов структуры «ведущего» объекта, а управляемое преобразование реализуется только в соот- 
ветствующих технологических процессах. 

2. Сопряжение «ведомого» и «ведущего» объектов взаимодействия наложением связей своего функци- 
онального назначения со стороны в общем случае группы элементов структуры «ведущего» объекта на функ- 
циональные элементы в структуре «ведомого» объекта позволяет определить метод преобразования «ведомого» 
объекта и свойства сопрягаемых объектов. 

3. Технологические операции информационной технологии решения задач стадий КПП и ТПП распре- 
деляем на основные и вспомогательные операции. Основные операции распределяем в два направления — тех- 
нологическое и конструкторское. В основных операциях рассматриваем виды обеспечения технологических си- 
стем организации рабочих машин и систем рабочих машин для обработки резанием. Во вспомогательных опе- 
рациях рассматриваем виды обеспечения технологических систем управления рабочими машинами для обра- 
ботки резанием. Вспомогательные операции конструкторского направления целесообразно закрепить за соот- 
ветствующими отделами стадии КПП. Виды обеспечения технологических систем организации и управления 
определяем в информационной и логической связи в составе: организационное, методическое, логическое, 
нормативное, информационное, программное, техническое. 
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4. Объектами преобразования в технической подготовке механообрабатывающего производства явля- 
ются материальные объекты {детали, сборочные узлы изделий } стадии КПП и {исходные заготовки, заготовки} 
стадии ТПП. Соответствующие материальным объектам стадий КПП и ТПП геометрические объекты являются 
идеальными. 

Структура объектов стадии КП рассматривается в соответствии с составом рабочих функций, испол- 
няемых каждой из четырех групп элементов структуры [15]. Рабочие функции элементов структуры в этой ра- 
боте определены объединением функций базирования с другими функциями в возможные их сочетания. Состав 
функций: базирование, направляющие, передача крутящего момента, деление и фиксация. 

Структура объектов {исходные заготовки, заготовки} стадии ТПП рассматривается в соответствии с 
составом рабочих функций, исполняемых двумя «ведомыми» объектами взаимодействия в соответствующих 
технологических схемах. Состав функций: базирование, базирование и геометрическое формообразование эле- 
мента конструкции. 

5. В основу познания геометрических объектов положено распределение основных элементов (инте- 
грации, дезинтеграции) структуры объектов производства (исходных заготовок, заготовок) соответственно и 
параметров их проектного качества по маршруту преобразований резанием на группы: 

» линейные и угловые размеры элементарных простых и налагаемых (простых, сложных) элементов 
формы; 

® макро геометрическая форма элементов формы; 

е микрогеометрическая форма (высота микронеровностей) элементов формы; 

е позиционные связи для исполнения функций взаимного расположения между центрами систем ко- 
ординат функционально различных групп элементов формы в структуре одного объекта и различных объектов, 
а также между координатными «геометрическими» осями «основных» и «налагаемых» простых и сложных 
элементов формы. 

» функционально различные отношения для исполнения функций взаимного расположения, например 
(совмещение, перпендикулярность, параллельность) между координатными «геометрическими» осями элемен- 
тов формы в различных их группах структуры одного объекта и между координатными «направляющими» ося- 
ми технологических схем взаимодействия различных объектов. 

6. Информационные модели систем организации рассматриваем на основе следующих принципов: 

» вне времени, в пространстве в общем случае совокупности информационно и логически связанных 
систем координат; 

е идеально в условиях статической настройки технологических схем взаимодействия информацион- 
ных объектов без применения каких-либо сил, но с учетом фактора времени в условиях динамической настрой- 
ки процессов управления преобразованиями соответствующих реальных «ведомых» объектов; 

е с учетом функционально различных погрешностей параметров проектного качества основных эле- 
ментов интеграции, дезинтеграции структуры реальных «ведомых» объектов. 

7. Технологию информационного логического моделирования решений задач практики методами и со- 
ответствующими средствами рассматриваем на основе следующих принципов. 

7.1. Оперируем с функционально определенной системой понятий, каждое из которых имеет свое фор- 
мализованное обозначение. Решение задач проводится на основе распределения структуры информационных 
моделей для основных базовых объектов знаний в две части информации: инвариантная, типовая объектно- 
ориентированная параметрическая. Информационные модели базовых объектов знаний более высокого уровня 
организованы на основе типовых объектно-ориентированных параметрических частей. 

7.2. В типовых объектно-ориентированных параметрических частях информационных моделей базовых 
объектов знаний всех типов решаются задачи безусловного обеспечения проектного качества основным эле- 
ментам (интеграции, дезинтеграции) структуры «ведомых» объектов преобразования. Решение задач всегда 
конкретно в пределах записи информации. Проводится в соответствии с интегральным алгоритмом в техноло- 
гии отображения и преобразования базы знаний предметной области в иерархии структуры классификации ее 
базовых объектов знаний по семи уровням структуры; в соответствии с графом структуры базового объекта 
знаний какого-либо типа. 

7.3. В базах данных решений определенных на основе инвариантных частей информационных моделей 
основных базовых объектов знаний в общем случае в пределах преобразований структуры, рабочих функций и 
параметров различных свойств решаются задачи оптимизации материальных и трудовых ресурсов. 

7.4. Используем средства системного анализа знаний: концептуальные понятия, теорию множеств, тео- 
рию графов. Концептуальные понятия — это «отображение», «преобразование», «структура», «множество», др. 
Графические средства — формализованные обозначения понятий встраиваемых в структуру базовых объектов 
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знаний. Теория множеств — раздел математики, включающий понятие «логический оператор» для высказыва- 
ния утверждений наложением связей между формализованными понятиями и т.д. Теория графов — раздел ма- 
тематики, изучающий, в частности, «графы структуры», используемый «диаграммы Венна» для наложения свя- 
зей между функционально едиными объектами рассмотрения. 

8. Информационное логическое моделирование базовых объектов знаний функционально различных 
типов какой-либо предметной области, в основу которого положены принципы системного анализа для нала- 
живания непрерывного, гибкого интегрального алгоритма в технологии (методы, соответствующие средства, 
алгоритмы, операции) отображения и преобразования знаний, выявляет все «огрехи» аморфного (вне структу- 
ры) отображения и преобразования знаний. Здесь важны целевые свойства понятия «логическое моделирова- 
ние» и функции знаний, следующие из его формулирования. 

9. Математические модели процессов управления рассматриваем на основе следующих принципов: 

е в соответствии с алгоритмом в технологии решения задач предметной области знаний; 

® с учетом информационных и логических связей органического единства систем организации и про- 
цессов управления, параметров времени и пространства; 

е не идеально в условиях динамической настройки технологических процессов взаимодействия ре- 
альных объектов технической, физической, химической, биологической природы с применением различных 
СИЛ. 

Определяемые экспериментальными измерениями различные погрешности параметров проектного ка- 
чества основных элементов (интеграции, дезинтеграции) структуры реальных «ведомых» объектов в процессе 
их преобразования подлежат учету с последующим нормативным их упорядочиванием до расчетных предель- 
ных допускаемых отклонений точности для использования в справочной документации в виде рекомендуемых 
технических условий на проектирование. 

10. Технологию математического моделирования решений задач практики рассматриваем на основе 
следующих принципов: 

» оперируем с функционально определенной системой параметров, каждый из которых обозначен 
своим символом; 

е с использованием различных методов и соответствующих средств математического аппарата; 

® решение задач всегда конкретно в пределах записи информации, проходит как преобразование па- 
раметров (параметрические модели). 

Результаты исследования. Под управлением функции целеполагания человека [16] понятия «про- 
странство — время», «системы организации — процессы управления», «логическое моделирование — матема- 
тическое моделирование» рассматриваем на инвариантной основе относительно каких-либо предметных обла- 
стей знаний в их системной информационной и логической связке. Идентификатором их связи, обладающим 
свойством органического единства, отрицающим раздельное рассмотрение, является понятие «содержание пре- 
образования «ведомых» объектов» в пределах определенных технологических элементов управления [17]. 

Технология системной инженерии, определенная автором для ТПП механообработки в практике дого- 
ворных работ как информационная технология познания базы ее знаний, может быть распространена на давле- 
ние, сварку, функционально различный менеджмент, компьютерную обработку информации. 

Обсуждение и заключение. Стадия ТПП механообработки, безусловно, крупное сложное информаци- 
онное изделие, включающее в себя как функционально различные системы организации составляющих ее ин- 
формационных изделий, так и процессы управления в них. Множественные попытки решения задач автомати- 
зации ТПП по всему циклу ее практики оказывались неудачными по причине логически неформализованного 
представления исходных данных и наличия существенных трудностей, характерных для крупных сложноорга- 
низованных автоматизированных систем, при реализации проектных решений с использованием только мате- 
матического аппарата. Исходя из логической незавершенности традиционных описательных знаний, переход к 
технологии их отображения и преобразования в знания информационной природы для достижения достаточно- 
го уровня автоматизации решения задач практики существенно затруднен или невозможен. Для налаживания 
непрерывного и гибкого интегрального алгоритма в детерминированной технологии решения задач высокого 
уровня автоматизации необходимы методы и инструментальные средства системной инженерии информацион- 
ного логического моделирования. 
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Моделирование и численный анализ влияния реакции диссоциации (рекомбинации) 
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Введение. Статья посвящена теоретическому изучению 
переноса ионов бинарной соли с учетом реакции диссоци- 
ации (рекомбинации) воды. Цели исследования: построе- 
ние математической модели; разработка алгоритма чис- 
ленного решения краевой задачи, соответствующей мате- 
матической модели; развитие теории подобия, включая 
переход к безразмерному виду с использованием харак- 
терных величин; определение физического смысла триви- 
альных критериев подобия; нахождение нетривиальных 
критериев подобия; построение и анализ вольтамперной 
характеристики (ВАХ). 

Материалы и методы. При теоретическом исследовании 
и численном анализе переноса ионов бинарной соли учи- 
тывается реакция диссоциации (рекомбинации) воды. При 
этом используются уравнение теплопроводности и мате- 
матическая модель электродиффузии одновременно четы- 
рех видов ионов (двух ионов соли, а также ионов Н\* и 
ОН”) в диффузионном слое электромембранных систем с 
идеально селективной мембраной. 

Для дифференциальных уравнений первого порядка ста- 
вится сингулярно возмущенная краевая задача. В уравне- 
нии для напряженности электрического поля правая часть 
не зависит от напряженности. При численном решении 
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ваются недостающие краевые условия для напряженности 
электрического поля. 

Результаты исследования. Разработаны новая математи- 
ческая модель, алгоритм численного решения краевой 
задачи, программное обеспечение. Проведен численный 
анализ и определены фундаментальные закономерности 
переноса ионов соли с учетом реакции диссоциации (ре- 
комбинации) молекул воды, температурных эффектов и 
Джоулева разогрева. Построена и проанализирована ВАХ. 
Обсуждение и заключение. Рассмотрен перенос ионов 
бинарной соли через диффузионный слой у катионооб- 
менной мембраны. Предложена математическая модель 
названного процесса, которая учитывает температурные 
эффекты, обусловленные реакциями диссоциации (реком- 
бинации) молекул воды и Джоулевым нагревом в раство- 
ре. Установлены основные закономерности переноса 
ионов соли с учетом реакции диссоциации (рекомбина- 
ции) молекул воды и температурных эффектов. Темпера- 
турные эффекты от реакции диссоциации (рекомбинации) 
и Джоулева разогрева в ОЭН практически незаметны (ис- 
ключение — область рекомбинации, ОР). Джоулев нагрев 
в области пространственного заряда (ОПЗ) на два порядка 
больше охлаждающего эффекта реакции диссоциации 
воды. При рекомбинации в ОР выделяется примерно 
столько же тепла, сколько при Джоулевом нагреве в рас- 
ширенной ОПЗ. Однако из-за малых размеров ОР влияние 
этого тепла незаметно. Значит, можно считать, что есть 
только один источник тепла на межфазной границе в ОПЗ, 
который благодаря своему заметному размеру обусловли- 
вает значительное повышение температуры во всем диф- 
фузионном слое. Отсюда следует, что возможно возник- 
новение и развитие гравитационной конвекции. Общие 
выводы, следующие из полученных результатов, открыва- 
ет возможность интенсификации процесса переноса ионов 
соли в электродиализных аппаратах. 


Ключевые слова: мембранные системы, ионообменная 
мембрана, пространственный заряд, расширенная область 
пространственного заряда, диссоциация (рекомбинация) 
молекул воды. 
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Введение. Интенсивный токовый режим предполагает использование токов, в несколько раз превыша- 


ющих предельный электродиффузионный. В этом случае в мембранных системах возникают вторичные (или 


сопряженные) явления концентрационной поляризации: 


— пространственный электрический заряд занимает макроскопическую область, сопоставимую с тол- 


щиной диффузионного слоя; 


— меняется рН раствора [1-3], что объясняется расщеплением воды; 


— формируются микроконвективные течения и пр. 


Учет влияния реакции диссоциации (рекомбинации) молекул воды важен для понимания процессов пе- 


реноса в электромембранных системах, поскольку ряд авторов считает, что появление новых носителей заряда 


Н* и ОН` может обусловить уменьшение или даже исчезновение пространственного заряда, являющегося ос- 


новой других механизмов переноса (например, электроконвекции). Реакция диссоциации молекул воды — эн- 
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дотермическая, а реакция рекомбинации — экзотермическая. Центры реакционных областей разделены в про- 
странстве. Это вызывает неравномерное распределение температуры и, следовательно, влияет на все физиче- 
ские свойства воды: коэффициент диссоциации, теплоемкость и пр. Кроме того, неравномерное распределение 
температуры может вызвать гравитационную конвекцию. Таким образом, исследование температурных эффек- 
тов, связанных с реакциями диссоциации и рекомбинации молекул воды, является актуальной задачей [4]. В 
данной работе ограничимся учетом влияния неравномерного распределения температуры на коэффициент дис- 
социации молекул воды. 

Материалы и методы 

1. Математическая модель 

1.1. Система уравнений 

Рассмотрим одномерный стационарный процесс переноса для однозарядных ионов соли в диффузион- 
ном слое Нернста. Учтем при этом реакции диссоциации (рекомбинации) воды и температурные эффекты, свя- 
занные с Джоулевым разогревом раствора. Такой случай описывается системой уравнений [5]: 


: —Е аф АС 
Л1 т =. (1) 
: Е аф ас 
]2 и (2) 
. —Е аф ас 
]з ао (3) 
: Е аф ас 
14 = ть 44 ах 4 =. (4) 
а?ф _ -Е 
ах? = =, (С1 — (2 + (3 — С), (5) 
ыы а 
-- =. +А; ‚Е=1....,4, (6) 
а?т 
к +Е-1 + В(СзС4 — Ки(Т)) = 0. (7) 
Здесь В. = В» = 0,В; = В. = КаСн»о — К„СзСа = К, (К, — СзС4), Ки(Т) = К ‚еР(Т-То), + = 1,2 — индексы ионов 
соли; { = Зи [= 4 — соответственно индексы ионов водорода Н* и гидроксила ОН”; Е = - — напряжен- 
ность электрического поля; С;, д, О; — соответственно концентрация, поток, коэффициент диффузии 1-го иона; 


Ка(Г) — константа скорости диссоциации воды; К.. — константа скорости рекомбинации ионов водорода Н* и 

гидроксила ОН”; К„— константа равновесия; =. — диэлектрическая проницаемость раствора; Е — число Фа- 

радея; Г — общая плотность тока; К — коэффициент теплопроводности; ро — плотность раствора; с, — удель- 

ная теплоемкость раствора; @ = Ё [Г — плотность источников тепла, связанного с прохождением электриче- 
Касн2о 

ского тока через раствор (Джоулев разогрев раствора); В = акК.., О = В (с; 4 — ие) = “К, 16: 5. (Г)) — 


плотность источников и стоков тепла, связанного с реакциями рекомбинации и диссоциации молекул воды. 
Общая плотность тока Г равна: [ = [$ + [,,. 
5 =Е (1 -]2), № =Е: Цз — 44). (8) 
Здесь [; — плотность парциального тока по ионам соли; [, — плотность парциального тока по ионам 
ВОДЫ. 
Из (6)-(8) следует: и = соп5$®, }› = соп5®, [ =ЁЕ` (1 -]2) = соп5$, [, = с0п5$6, но }з = с0П$Е, 4 = 
К 
СОП5Е, а? = ти коэффициент температуропроводности раствора. 
обр 


1.2. Краевые условия 

Будем считать, что х = 0 соответствует глубине раствора, где выполняется условие электронейтраль- 
ности С: — С. + С. — С, = 0, ах = 6 соответствует условной межфазной границе «раствор — мембрана». Та- 
ким образом, 5 — толщина диффузионного слоя Нернста. 

В этой работе для определенности рассматривается катионообменная мембрана (КОМ), которая для 


<“ „(14 
простоты будет предполагаться идеально селективной . = 0). 





Рис. 1. Схема диффузионного слоя у катионообменной мембраны 


Чубырь Н. О. и др. Моделирование и численный анализ влияния реакции диссоциации (рекомбинации) молекул воды 


Сйируг М. О., её. 1. Мо4аейп? апа питетса! апа[уя15 ор йе ерресй ор @55осаноп/тесотЬтапоп оГижиег тоесшех оп ше гапурот ор зай 





1) Краевые условия в глубине раствора (х = 0) 
Положим, С1(0) = Сло, С2(0) = Соо, Сз(0) = Сзо, С4(0) = —— "> причем Сло + Сзо = С2о, Т(6) = То. 


2) Краевые условия на межфазной границе И (х=б) 

Полагаем, С, (6) = С\т, С3(6) = Сзт. При этом сумма Ст + Сът сопоставима с табличным значением 
обменной емкости катионообменной мембраны. Плотность потока (5) = мт ионов ОН `` определяется ката- 
литической реакцией диссоциации молекул воды в тонком поверхностном слое катионообменной мембраны. 
Для температуры используется условие баланса потоков в точке х = 6 [1]: 

НКТ = 6(5)6 + 0(6)5 +49 

2. Переход к безразмерному ви 

2.1. Характерные величины и формулы перехода к безразмерному виду 

В качестве характерной концентрации Су возьмем концентрацию катионов соли С10о. Реальные значения 
С1о находятся в диапазоне 10-?М -— 1М. В качестве характерной температуры возьмем комнатную температуру 
То = 298К. В качестве характерной длины рассмотрим толщину диффузионного слоя 6, который имеет порядок 








10 %м. Обозначим Ду = о — коэффициент диффузии электролита (МаС), Фо = = — тепловой потенци- 
1 2. 
ал [5]. 
о к безразмерным величинам: 
хо = И фе № Е = —6 ев, 19) = — 2% = > 1% = —— Е =”,  ^ — _- _ — о = 
0 Фо ах пр пр” тр Сто Р:С1о 
ЛЕ — (и) — — бб - 
‚= 1,2,3А, т =- 1, = 
Е т 1 пр 5 
(м) —_ КиСГ) _ (то — —_ Км —_ 
Тогда Ки” (Г) = и = уе (т 1). где у@2 2 р = ЬТь. 
О АИ ты 
Положим, & С оЕ252° т 


2.2. Переход к безразмерному виду в уравнениях 
Ориентируясь на [6], перейдем к безразмерному виду в системе уравнений (1)—(6): 


(и) 
аб 1 _ (м) .(и) 
= С: Е д“, (9) 
(и) 
а(> ^ _ _ (и) (и) _ и 
о (10) 
0 
асз— (и) (и) _ (4) 
и Е (11) 
(и) 
ао _ _ „(м . (и) 
пе о. (13) 
0 р@д ры ыы (Рева) —_ 6. (14) 
а _ (м) _ „ад, 0) _ пы) 
В о (15) 


Запишем общую плотность тока в виде 1%) + 1 = 1%. 


Здесь 


16) № вы в” 


(1) 


= 


= 0 и, следовательно, 
1%) — 2“ (и) 


причем для идеально селективной мембраны ]> 








Л 
При переходе к безразмерному виду уравнение теплопроводности (7) приобретает вид: 
р (и) (ти) — 
ре ЕО 199 + ВО (СС — удев" 1-9) = 0. (16) 
Здесь К(Э = 60 Куб” С1о 
РоС1оВ РовТо ° 


2.3. Физический смысл и оценка величин тривиальных критериев подобия 

Параметры =), В р а, ко, К 2, В05, Ь@®) являются тривиальными критериями подобия. Опре- 
делим их физический смысл и значения при характерных значениях входных параметров задачи. 

1. 5% = 15. 

2. Параметр у представляет собой квадрат отношения концентрации ионов Н* или ОН” в нейтраль- 
ном растворе (ФН = 7) к концентрации электролита. Его значения меняются в пределах от 10—10 до 10-18, ион 
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можег считаться малым параметром. 





2 

ВТо= 1 ВТЕ 

3. = = о. =” (4) ‚ где [1 = | > — дебаевская длина. Таким образом, параметр =@9 является 
10 10 


удвоенным квадратом отношения дебаевской длины [1 к ширине диффузионного слоя б. Значения = меняют- 
ся в пределах от 10-11 до 1077, и этот параметр можно считать малым. 





2 

ВТоЕу К ВТо=у" 2КуС [ р 

Е бы Сы (№) ‚ Здесь [.. = `_. Данная величина пропорциональная ширине 
Е? Ро 2б10Е? Ро Е бо 


области рекомбинации (ОР). Таким образом, параметр а@ — это квадрат отношения дебаевской длины к ши- 


рине ОР [., . Параметр а = 14 ив рамках данной задачи является универсальной постоянной, поскольку не 


[ 
зависит от входящих параметров задачи. Таким образом: _ = \/а 9 = 3,74, т.е. ширина ОР составляет при- 
т 


мерно 37 % от ширины квазиравновесной области пространственного заряда (ОПЗ, дебаевская длина). 
2.4. Нетривиальные критерии подобия, физический смысл и оценка величин 
В уравнении (16) безразмерные коэффициенты К@9 и В© являются тривиальными критериями подо- 
бия. Любое соотношение между тривиальными критериями называется нетривиальным критерием подобия. 
Ниже выразим К и В через другие известные критерии подобия. 
К а’роср _ а? Роср 2 


Рос © 
1) Рассмотрим КО = = ОР. ОР =. ОР Здесь — — число Льюиса — Семенова. Введем 
Роб1оК  Боб10о По С1оВ СтоК Бо 





Со 1 
в рассмотрение безразмерное число и = ЕЕ ‚ Тогда К@9 =. -. 
оСр и 


Число и обозначает относительную удельную объемную теплоемкость, изменяется в пределах от 10-8 
до 10“, и его можно рассматривать как малый параметр. 








ак!-52 С а5? 2К..С а 52 52 а 
2) Рассмотрим параметр В = = "^_1® = —5— = Атп--. Здесь Ати = — = 22,8 — 
Ровто — 2ЕТо Бо 2ЕТо о 212 ВТо 
те10 


безразмерная величина, являющаяся тепловым аналогом критерия Аррениуса. Она характеризует чувствитель- 


2 
> 1 
ность количества тепла протекания химическои реакции к изменению температуры. Так как (>) — 


нада 2 ово ди 


Таким образом, уравнение теплопроводности (16) имеет следующий безразмерный вид: 








ат 
&Ёе 
3) Вычисления показывают, что 
= К» В _ КиЕЕ? 6" В _ Ку 62 
Сто РоСбр ВТо=у РоСр РоСрТоЕт 


Ки„Е262 
РоСрТо т. 


терий подобия, который показывает соотношение 


2 
1 и е 
Поскольку #09 = 2 (4) ‚ то М = = 2,5 . 10-3. Следовательно, получаем нетривиальный кри- 


и. 


ун = Па. 
Аналогично 
ЕН = Пал 
Здесь 
ВАТ —14 
ПЕ 30.10 17 
а,1 росрЕ? 52 , ( ) 
— 2 
= Ма 2Е > (18) 
2464 К 
М2 = — 
В2Т2 Е? 


Заметим, что п2 = Па1 Ма. 

2.5. Безразмерный вид системы уравнений 

Таким образом, безразмерная система уравнений, описывающая модель, имеет вид (индекс и опущен 
для простоты записи): 


ас ‚ ас ‚ ас у 

к = СЁ — Лох Я Е — рта, 2 СЁ — Л 
а. в ас ; 

ЕБз--)з = а(уе?Т-® — С.С.), = р 


а. — Е 
Е. а = а(уе?-Ъ — С.С,), 52. = С, — С> + С3 — Са. 


а?т — 
Ге > + ЕНЕТ + цаАтт (СзС. — уе 75) = 0. 
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Система уравнений содержит три малых параметра &, ц и у, связанных соотношениями (17) и (18). Па- 
раметр = входит в уравнения сингулярно, т. е. как коэффициент при производных, а параметр у — регулярно, 
поэтому система уравнений является одновременно сингулярно и регулярно возмущенной. 

2.6. Переход к безразмерному виду в краевых условиях 

Безразмерные краевые условия имеют вид: 


Су|х=о = 1, С>|х=0 = 1 Сзо, Сз|х=о = Сзо, Са|х=о = Со = Г(0) = 1, 
С. х=1 = Слть Сз[=а = Сзт» Ла|х=1 = Ллть 
‚ а Я. 
Ка ГС) = НО4Ба чел. 
Замечание 1. В силу непрерывности справедливо: 
Лу = Эз]зт — Э4Лат = Оз]зо — Вало. 
3. Алгоритм численного решения 
3.1. Переход к краевой задаче для системы уравнений второго порядка 
Из-за наличия малого параметра = в уравнениях (12), (14) и (15) краевая задача (9)—(16) относится к 
сингулярно возмущенным, поэтому для численного решения система неудобна [7]. Кроме того, дифференци- 
альные уравнения имеют первый порядок, а для них ставится краевая задача. В уравнении (15) для напряжен- 
ности электрического поля правая часть не зависит от Е. Этим обусловлены проблемы с устойчивостью при 
численном решении дискретизированной системы уравнений методом Ньютона — Канторовича [8]. Поэтому 


краевую задачу (9)—(16) следует привести к удобному для численного решения виду. Для этого перейдем к си- 
стеме уравнений второго порядка: 








аб: .. @. а2с› _ -а а? _ а — а(уе (Т--с3с.) 
ах? ах (С.Е), ах? ах (С›Е), ах? ах (СзЕ) =Рз ? 
ас, _ а — ат №364) „аЕ _ пр 96 _ 964 
ах? ах (С4Е) =Ра | 2 — Е(С, ы С») тт ах ах 
а?т Ь(Т-1)\\ — 
Ее + НЕТ +: цаАти (С3 С. — уе ) = 0. (19) 


В связи с переходом к системам уравнения второго порядка необходимо дополнить исходные краевые 
условия недостающими краевыми условиями для функции Е. Из краевых условий следует, что в точке х = 0 
выполняется условие электронейтральности С1о — (50 + Сзо — Со = 0. В силу непрерывности с некоторой точ- 
ностью в окрестности точки х = 0 выполняется С; — С. + С. — С, = 0, поэтому Еу может быть вычислено из 
условия электронейтральности. Складывая уравнения для (:, 1 = 1, ...,4 и используя условие электронейтраль- 


+4 (С1(0) — С2(0)). 


_ = (С, — С> + С; — (4) = 1. 


х= 





ности (10 — С50 + Сзо — Со = 0, можно получить соотношение ЁЕ|„-о = 


АЕ 
Из уравнения для Ё прих = 1 имеем: Е 
Хх 


Результаты исследования 
3.2. Результаты численного решения 
Для анализа результатов численного решения наряду с графиками искомых функций С}, Е, }з, ]л, Г рас- 
смотрим графики функций р(х) = С, — С, + С; — С4, р(х) = С.С, — уе Ч-Ъ. Функция р(х) характеризует рас- 
пределение плотности зарядов, а р(х) — отклонение от равновесия реакции диссоциации (рекомбинации). 
Глпе ОтарН:С3*С4-отатта*ехр(Ь*(Т-1)) 


х10 
1.6 
1.5 
1.4 
1.3 
1.2 
1.1 


0.9 
0,8 
0.7 
0.6 
0,5 
0.4 
0,3 
0.2 
0.1 


С3*С4-ватта*ехр(Ь*(Т-1)) 


0 0,1 0,2 0,3 0.4 0,5 0,6 0,7 0,8 0,9 1 
Агтс еп 


а) 
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С1-С2-+С3-С4 


Переп4еп{ уапае Т 


С3*С4-ватта*ехр(Ь*(Т-1)) 
= 





Глпе ОгарБ:С1-С2+С3-С4 


0.006 
0.005 
0.004 
0.003 
0.002 


0.001 


0.001 
—0.002 
0.003 


—0.004 
0 0.1 0.2 0.3 0,4 — 0,5 0.6 0.7 0.8 0,9 
Атс еп 


6) 


Глпе ОтарН:С3*С4-отатта*ехр(Ь*(Т-1)) 


х10- 


—1,2 
—1,4 
—1,6 
—1,8 


0.5 0.81 0,82 0,83 0.84 0,85 0.86 0,87 0,88 0,89 


Агс еп 


в) 
Глпе ОгарВ: Оерепдет{ уапае Т 


1.0028 
1.0026 
1.0024 
1.0022 

1.002 
1.0018 
1.0016 
1.0014 
1.0012 

1.001 
1.0008 
1,0006 
1.0004 
1.0002 





0 0.1 0.2 0.3 0.4 0,5 0.6 0,7 0,8 0,9 1 
Атс еп 


г) 


Рис. 2. График функции р(х) = С, — С + С; - С, (10: 0,93] (а); график функции р(х) = СзС. — уе Т-1(6); 
график функции р(х) = СС. — уе -Ъ в ОПЗ (в); график температуры (2) 


Из рис. 2 а следует, что диффузионный слой делится на три области. Первая — область электроней- 


тральности (ОЭН), где р(х) = 0. Вторая — область пространственного заряда (р(х) > 0). Третья — промежу- 
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точный слой между ними, где плотность распределения зарядов имеет ярко выраженный локальный максимум. 
Сама ОПЗ делится на две части: расширенную и квазиравновесную (она очень мала и примыкает к межфазной 
границе). В квазиравновесной части ОПЗ распределение плотности зарядов подчиняется экспоненциальному 
закону (на рис. 2 не приведено). 

Общий вид графика функции р(х) (рис. 2 6) показывает, что отклонение от равновесия происходит: 

— в области рекомбинации, которая находится внутри ОЭН; 
— в области пространственного заряда (рис. 2 в). 

В первом случае р(х) > 0, то есть рекомбинация преобладает над диссоциацией С3С. » уе? ("-2. Во 
втором случае, наоборот, р(х) < 0, и диссоциация преобладает над рекомбинацией в ОПЗ, а также в промежу- 
точном слое. Температура (рис. 2 г) увеличивается линейно на 0,0028 в безразмерном виде (или 0,0028 : 298 = 
0,83 градуса). Увеличение зависит от общей плотности тока или скачка потенциала в диффузионном слое и 


может составить до десятков градусов. Очевидно, что все члены уравнения (19) по сравнению с первым малы, 
2 


а“т В: 
то есть вторая производная ав = 0. Из этого следует, что температура является линеинои и определяется крае- 
Хх 


выми условиями. 
Глпе ОгарВ: Оереп4деп{ уамае С1 Глпе ОтарВ: Вереп4еп{ уапа е С2 


Глпе ОгарВ: Оереп4деп{ уамае СЗ Глпе ОтарВ: Вереп4еп: уапа ]е С4 


—— С1 (Ма+) 
СС) 
—— С3(Н+) 
—— 64 (ОН-) 
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Рис. 3. Общий вид графиков функций Ст, (>, С5, С (а); графики функций С\ С в окрестности нуля (6); графики функций С3, 
С. в окрестности нуля (в); график функций Ст, С, С3, С. в области пространственного заряда (г) 


Как следует из рис. 3 а, концентрация натрия линейно падает в ОЭН и становится малой в ОПЗ, но она 
выше, чем концентрация остальных ионов (рис. 3 г). Центром ОР является точка пересечения графиков концен- 
трации Ма" и СГ и соответственно ОНТ и Н*. Правее центра концентрация С[` убывает быстрее, чем концен- 
трация Ма". Возникающий из-за этого дефицит отрицательно заряженных ионов компенсируется быстрым ро- 
стом ионов гидроксила, и сохраняется условие электронейтральности (рис. 3 а). 
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Рис. 4. Общий вид графика функции Ё за исключением межфазной границы (а); график функции Ё вблизи ОР (6); 
график функции Ё вблизи межфазной границы (в), вольтамперная характеристика (г) 
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Из данных рис. 4 следует, что в ОЭН напряженность электрического поля быстро увеличивается в де- 
сятки раз (рис. 4 а). Особенно высокая скорость роста напряженности электрического поля отмечается в ОПЗ 


1 2 < 
при достижении величин порядка --. В ОР имеется внутренний погранслой по напряженности электрического 
Е 


поля Ё (рис. 4 6). Это связано с тем, что меняются носители заряда: до ОР заряд определяется ионами Ма*, СГ, 
Н*, а после — ионами Ма", СГ, ОН". Появление «ступеньки» связано с тем, что коэффициент диффузии у 
водорода примерно в два раза больше коэффициента диффузии гидроксила. ВАХ неограниченно растет с уве- 
личением скачка потенциала в запредельном режиме [; > 2 (рис. 4 г). Этот рост происходит значительно мед- 
леннее, чем экспериментально наблюдаемый — значит, наряду с реакцией диссоциации (рекомбинации) воды 
существует и другой механизм сверхпредельного переноса. Общепринято, что таким механизмом является 
электроконвекция [9-11]. В связи с этим возникает проблема оценки и взаимодействия реакции диссоциации 
(рекомбинации) и электроконвекции. 
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Рис. 5. Графики потоков ионов водорода ]з (а) и гидроксила ]^ (6) 

Из данных рис. 5 видно, что поток ионов водорода положительный. Он постоянен слева от ОР и прак- 
тически равен нулю справа от ОР вплоть до ОПЗ. В ОПЗ поток ионов водорода линейно увеличивается в связи 
с интенсивной диссоциацией воды. Поток ионов гидроксила левее ОР практически равен нулю, а правее — от- 
рицателен и постоянен до ОПЗ. В ОПЗ поток ионов гидроксила так же, как поток ионов водорода, линейно рас- 
тет. В квазиравновесной ОПЗ поток ионов гидроксила резко меняется и удовлетворяет граничному условию 
(1) = „(т). Область рекомбинации является внутренним погранслоем для потоков водорода и гидроксила — 
здесь они резко меняются. 

Обсуждение и заключения. В данной работе рассмотрен перенос ионов бинарной соли через диффу- 
зионный слой у катионообменной мембраны. Предложена математическая модель названного процесса, кото- 
рая учитывает температурные эффекты, обусловленные реакциями диссоциации (рекомбинации) молекул воды 
и Джоулевым нагревом в растворе. Эта модель состоит из 9 обыкновенных дифференциальных уравнений с 
соответствующими краевыми условиями. Разработан алгоритм численного решения краевой задачи, приведено 


ее численное исследование. Установлены основные закономерности переноса ионов соли с учетом реакции 


Чубырь Н. О. и др. Моделирование и численный анализ влияния реакции диссоциации (рекомбинации) молекул воды 


Сйируг М. О., её. 1. Мо4аейп? апа питетса! апа[уя15 ор йе ерресй ор @55осаноп/тесотЬтапоп оГижиег тоесшех оп ше гапурот ор зай 





диссоциации (рекомбинации) молекул воды и температурных эффектов. 

Показано, что ВАХ неограниченно растет с увеличением скачка потенциала в запредельном режиме, 
причем значительно медленнее, чем в экспериментах. Следовательно, наряду с реакцией диссоциации (реком- 
бинации) воды существуют и другие механизмы сверхпредельного переноса (например, электроконвекция и 
гравитационная конвекция). Поэтому возникает проблема оценки взаимодействия реакций диссоциации (ре- 
комбинации), электроконвекции и гравитационной конвекции. 

Температурные эффекты от реакции диссоциации (рекомбинации) и Джоулева разогрева в ОЭН прак- 
тически незаметны (исключение — ОР). Джоулев нагрев в ОПЗ на два порядка больше охлаждающего эффекта 
реакции диссоциации воды. При рекомбинации в ОР выделяется примерно столько же тепла, сколько при Джо- 
улевом нагреве в расширенной ОПЗ. Однако из-за малых размеров ОР влияние этого тепла незаметно. Значит, 
можно считать, что есть только один источник тепла на межфазной границе в ОПЗ, который благодаря своему 
заметному размеру обусловливает значительное повышение температуры во всем диффузионном слое. Отсюда 
следует, что возможно возникновение и развитие гравитационной конвекции. 

В работе рассматривается раствор хлористого натрия, однако проведенное исследование справедливо 
для раствора любой бинарной соли. 

Общие выводы, следующие из полученных результатов, открывает возможность интенсификации про- 


цесса переноса ионов соли в электродиализных аппаратах. 
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Проблема математического конечноэлементного моделирования неоднородных 


деформируемых твердых тел с применением сканирования` 


В. Л.Зыонг”" 


' Иркутский национальный исследовательский технический университет, г. Иркутск, Российская Федерация 
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Введение. При математическом конечноэлементном моде- 
лировании используется усредненное значение механиче- 
ских характеристик материала деформируемых твердых 
тел. В авиа-, машиностроении, строительстве, медицине и 
других областях все шире применяются полимерные ком- 
позитные материалы и материалы природного происхож- 
дения. В последнем случае реальное изменение механиче- 
ских характеристик значительно отличается от усреднен- 
ного, следовательно, при использовании усредненных 
параметров для построения и анализа конечноэлементных 
моделей результаты могут существенно искажаться. 

В данной статье описано создание математических мето- 
дов исследования изменения механических характеристик 
материала неоднородных деформируемых твердых тел. 
Полученные таким образом результаты применены для 
построения конечноэлементных моделей и анализа их 
напряженно-деформированного состояния. 

Материалы и методы. В качестве неоднородных дефор- 
мируемых твердых тел рассмотрены материалы природ- 
ного происхождения и композиты. Для исследования из- 
менения механических характеристик материала разрабо- 
тан способ, основанный на использовании двух составля- 
ющих: пиксельной характеристики растровых изображе- 
ний сканирования компьютерным томографом и экспери- 
ментальных данных натурных испытаний стандартных 
образцов. 

Результаты исследования. Создан комплекс математиче- 
ских методов моделирования интерпретации растровых 
изображений сканирования компьютерным томографом, 
позволяющий проводить исследование любых сложных 
структур реальных деформируемых твердых тел. Резуль- 
таты используются при построении конечноэлементных 
моделей таких тел с учетом неоднородности механиче- 
ских характеристик материала. 

Анализ напряженно-деформированного состояния конеч- 
ноэлементных моделей тестовых образцов доказал точ- 
ность и сходимость численного решения метода конечных 
элементов при моделировании свойства неоднородности 
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механических характеристик материала. 

Обсуждение и заключения. Разработанный подход может 
быть применен для любых физических принципов скани- 
рования (рентгеновский, ультразвуковой, лазерный и др.) 
и для любых типов материалов, если информация, полу- 
ченная в результате сканирования, сформирована в виде 
цифрового (растрового) изображения. 


Ключевые слова: метод конечных элементов, деформи- 
руемое твердое тело, неоднородность, механические ха- 
рактеристики материала. 
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матического конечноэлементного моделирования неодно- 
родных деформируемых твердых тел с применением ска- 
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Введение. При конечноэлементном (КЭ) моделировании механические характеристики материала де- 
формируемых твердых тел (ДТТ) задаются в виде усредненного значения, используемого для всей модели. 
Например, механические характеристики стали определяются при испытании на растяжение и сжатие стан- 
дартных образцов правильной формы. Такой подход приемлем для расчетно-инженерной практики. Однако в 
современном авиа-, машиностроении, строительстве, медицине и других сферах все шире применяются поли- 
мерные композиты и материалы природного происхождения. Реальное изменение механических характеристик 
таких материалов сильно отличается от усредненного значения, следовательно, при использовании усреднен- 
ных параметров для построения и анализа конечноэлементных моделей результаты могут существенно иска- 
жаться. 

Результаты исследования 

1. Основная математическая зависимость моделирования неоднородности реальных деформиру- 
емых твердых тел. Для решения представленной проблемы в работе предлагается способ определения реаль- 
ного изменения механических характеристик материала, основанный на использовании: 

— результатов сканирования ДТТ; 
— данных натурных испытаний стандартных образцов [1]. 

Такой способ дает возможность распознавать любую сложную неоднородную структуру материала и 
применить эти данные для повышения точности и реалистичности при математическом конечноэлементном 
моделировании реальных ДТТ. 

На рис. 1 а представлен пример растрового изображения сканирования материала природного проис- 
хождения (человеческой кости). В ее сечении видно изменение механических характеристик в зависимости от 
значений индексов цвета (пикселей) [2]. На рис. 1 6 представлен результат сканирования композитного матери- 
ала в виде растрового изображения сечения образца правильной формы (сканирование выполнено 
В. Г. Толстиковым). Обозначены реальные изменения механических характеристик волокон и матрицы (клея) в 
структуре. Кроме того, сканирование показало возможность оценить качество композитных материалов (распо- 
ложение волокон, толщина слоя клея, внутренние дефекты и др.). 
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Рис. 1. Анализ растрового изображения сканирования материала природного происхождения 
(а) и композитного материала (6) 


Общая схема математического моделирования и алгоритм предлагаемой технологии сканирования 
представлена на рис. 2 |2]. 
1 


Деформируемое твердое тело 
с неоднородными свойствами 


2 Проведение сканирования 
компьютерным томографом 
3 Формирование пакета 
растровых изображений сечений 
4 Построение матрицы индексов цвета 
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Рис. 2. Математическое и конечноэлементное моделирование реального деформируемого твердого тела 





Анализ НДС КЭ модели 


Схема отражает полный комплекс модулей: от сканирования ДТТ до получения его КЭ модели и ана- 
лиза напряженно-деформированного состояния (НДС). При необходимости внутри каждого блока алгоритма 
обеспечивается цикличность операций. 

Сканирование проводится с использованием компьютерного томографа (КТ). При работе с деформиру- 
емыми объектами он позволяет формировать пакет снимков, каждый из которых представляет собой растровое 
изображение в определенном сечении [3, 4]. Стандартные форматы растрового изображения: ]рее, Ч1сот, Бтр, 
рие и др. Независимо от используемого формата растровое изображение необходимо представить как числовую 
матрицу для дальнейшей обработки. 

Чтобы определить зависимость механических характеристик (например, в виде модуля упругости) от 
значений индексов цвета, применяется специальный весовой коэффициент [1], устанавливающий зависимость 
между средним значением индексов цвета ие, и опытным (экспериментальным) значением Ёоп. Среднее значе- 
ние модуля упругости (Еоп), полученное в эксперименте на растяжение (сжатие) стандартных образцов, ис- 
пользуется при интерпретации механических изменений по пиксельной характеристике. 

Тогда весовой коэффициент, определяющий переход от индекса цвета пикселя к модулю упругости, 
определяется выражением вида: 
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Здесь и. — среднее значение индексов цвета пакета растровых изображений сканирования: 


У льл 
В (2) 


Однако для большого объема данных сканирования КТ ДТТ, имеющих высокую степень неоднородно- 
сти структуры, при определении среднего значения индекса цвета необходимо использовать зависимость мате- 
матического ожидания. То есть, и.„ определяется как среднестатистическая величина: 


пер = Е(п)= У трит), (3) 
1=1 
где и; — значение индекса цвета на интервале Г; р;(п;) — вероятность появления п; (значения индекса цвета) на 
интервале /. 

Пикселям, каждый из которых имеет соответствующий индекс цвета, присваиваются значения модуля 
упругости материала Е’. Для этого используется весовой коэффициент Ах: 

Е; =И;КК. (4) 
Здесь Ё; — значение модуля упругости, соответствующее значению номера индекса цвета и;; Ке — весовой ко- 
эффициент модуля упругости. 

Реализация зависимости (4) для КЭ модели выполняется решением следующих задач: 

— определение зависимости механических характеристик от значений индексов цвета; 
— определение механических характеристик материала в узлах и в конечных элементах КЭ модели. 

Изменение модуля упругости определяется на основе линейной или нелинейной зависимости от значе- 
ний индексов цвета. Для материала, имеющего невысокую степень неоднородности, зависимость может быть 
построена по линейному закону (4) с использованием единственного постоянного значения весового коэффи- 
циента модуля упругости. 

Однако в природе и технике такие материалы встречаются редко. Чаще всего в структуре материала 
ДТТ можно выделить несколько областей, каждая из которых имеет отличный от других диапазон изменения 
модуля упругости. 

Для таких материалов предлагается применять нелинейный закон, позволяющий использовать в каж- 
дой области структуры ДТТ свое значение весового коэффициента. В этом случае зависимость изменения мо- 
дуля упругости ДТТ от индексов цвета представляется в виде кусочно-линейной функции (рис. 3 а). 
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Рис. 3. Кусочно-линейная (а) и непрерывная (6) функция зависимости изменения модуля упругости ДТТ 
от значений индексов цвета 

Использование кусочно-линейной функции затруднительно с точки зрения эффективности вычисли- 
тельного алгоритма, в частности, по следующим причинам: 
— чем больше количество областей структуры ДТТ, тем больше количество весовых коэффициентов и тем 
сложнее вычисление модуля упругости для узлов и конечных элементов; 
—щ переходные значения индексов цвета между областями ДТТ или объектами представляют собой не одно 
значение, а диапазоны значений индексов цвета на растровом изображении (рис. 3 6). 

Для решения этой проблемы на основе кусочно-линейной функции можно построить нелинейную зави- 


симость в виде сплайна, например кубического [5]. Такая нелинейная зависимость модуля упругости ДТТ от 
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индексов цвета является универсальной при объединении нескольких зон с различным порядком механических 
характеристик (рис. 3 6). 

Далее. Интерпретация данных сканирования связана с интерполяцией и передачей значений индексов 
цвета от пикселей к узлам КЭ сетки. В самом сложном и общем случае узлы КЭ сетки лежат между соседними 
плоскостями с растровым изображением. Они существенно различаются по форме и размерам. Кроме того, в 
них появляются дополнительные внутренние контуры и другие геометрические особенности. В этом случае для 
интерполяции используется ряд математических зависимостей, определяемых относительно прямой, наклонной 
к оси / [6]. 

В [2] представлено более подобное описание математической зависимости и блок-схемы алгоритма мо- 
дуля определения механических характеристик материала в узлах и конечных элементах КЭ модели ДТТ 
(рис. 2, блок 7). 

Для использования в КЭ модели результаты определения неоднородности механических характеристик 
материала преобразуются в файл с расширением *.рс[. Для этого задействуют две функции языка Раёгап Сот- 
тапа Гапоцазе: «табета].стеже» и «@ететргорз _стеже» [7, 8]. 

2. Исследование точности и сходимости результатов анализа КЭ моделей с учетом моделирова- 
ния неоднородности механических характеристик. В рамках данной работы определены точность и сходи- 
мость численного решения методом конечных элементов (МКЭ) с применением представленной технологии. С 
этой целью использованы расчеты НДС костной ткани человека (фрагмент бедренной кости). 

Выбор ДТТ и его анализ непринципиальны, но обусловлены важными факторами: неоднородностью 
материала кости, индивидуальностью ее геометрии, а также высоким уровнем технологии и качества сканиро- 
вания КТ в медицине. Костная ткань хорошо изучена на практике, что позволяет обоснованно судить о кор- 
ректности ее математического моделирования. 

В первом исследовании используются образцы в форме параллелепипедов, вырезанные из фрагмента 
бедренной кости человека (рис. 4 а). Результаты исследования представлены на рис. 4 6 в виде доказательства 


точности и сходимости численного решения МКЭ относительно данных натурного эксперимента [9, 10]. 
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Рис. 4. Визуализация образцов: размеры и ориентация реальных образцов (а); график сходимости нормального напряжения 
в центре образцов (6), в продольном (1), окружном (2) и радиальном (3) направлениях 


Результаты показывают, что для достижения необходимой точности численного решения МКЭ 
требуется КЭ сетка плотностью три и более конечных элемента на 1 мм [11]. Очевидно также, что свойство 
неоднородности механических характеристик материала ДТТ может быть отражено набором конечных 
элементов с собственными модулями упругости и изотропной структурой материала. 





1 1,5 5 32 
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Рис. 5. Результат анализа НДС с применением КЭ типа НЕХ (а) и ТЕТ (6); график сходимости результатов (в) 
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Для изучения применимости двух различных типов объемных КЭ (гексаэдр и тетраэдр) проведено 
дополнительное исследование с использованием несовместных функций формы (рис. 5) [12]. Результаты 
показывают, что КЭ типа гексаэдр имеет лучшие показатели по точности и сходимости результатов, а также по 
ресурсным затратам. 

3. Моделирование реальных деформируемых твердых тел. Представленная технология применена к 
реальным объектам. С этой целью построено и проанализировано НДС КЭ моделей зубов с дефектом и 
пломбой из композитного материала. В одной модели представлен дефект в виде кариозной области (впадины) 
в верхней части эмали. Во втором варианте КЭ модели эта область заполнена пломбой из композитного 
материала. Обе модели построены с учетом неоднородности механических характеристик материала и 
индивидуальной геометрии. 

Каркасная модель зуба человека, представленная на рис.ба, построена по результатам 
сканирования [13, 14]. Изменение механических характеристик материала в объемной КЭ модели (рис. 6 6) 
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смоделировано в специальной программе [ 13, 14]. Входные данные: результаты сканирования и натурных 
испытаний. 
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Рис. 6. Построение и анализ КЭ моделей зуба человека с учетом неоднородности механических характеристик материала: 
каркасная модель (а); изменение модуля упругости зуба в объемной КЭ модели (6); расчетная схема (в); твердотельная 
геометрическая модель (г); генерация КЭ сетки (д); результат анализа НДС КЭ модели зуба с моделированием дефекта (е); 
результат анализа НДС КЭ модели зуба с моделированием композитной пломбы (ж) 


Расчетная схема представлена на рис. 6 в. На верхнем зубе моделируются дефект и композитная 
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пломба (рис. бгид). 

Результат анализа КЭ модели с дефектом (рис. бе) показывает, что кариозная часть зуба не 
воспринимает нагрузку. То есть на здоровую часть зуба приходится дополнительная нагрузка, и вероятность 
разрушения увеличивается. Результат анализа второй КЭ модели, после заполнения области кариеса пломбой 
из композитного материала (рис. 6 ж), показывает, что зуб практически полностью восстанавливает свои 
функциональные способности, так как картина напряжений идентична картине, полученной для здорового 
зуба [2]. 

Из рис. 6 ж видно, что максимумы напряжений наблюдаются в точках контакта на вершинах зубных 
бугров. Более высокая нагрузка воспринимается областью эмали и дентина. Внутренняя часть зуба менее 
нагружена. Данный результат сравнивается с расчетом [2]. При этом используется усредненное значение меха- 
нических характеристик эмали и дентина. Неоднородность структуры материала зуба не учитывается. Сравне- 
ние показывает более равномерное распределение напряжений в теле зуба при моделировании неоднородности 
распределения его механических характеристик. 

Обсуждение и заключения. Результаты исследования подтверждают точность и надежность представ- 
ленного метода КЭ моделирования на основе сканирования реальных ДТТ с учетом неоднородности механиче- 
ских характеристик материала. 

Представленный комплекс математических методов моделирования интерпретации растровых изобра- 
жений сканирования КТ позволяет проводить исследование любых сложных структур реальных ДТТ. Результа- 
ты такого моделирования используются при построении КЭ моделей с учетом неоднородности механических 
характеристик материала. 

Разработанная технология математического моделирования может быть применена для любых физиче- 
ских принципов сканирования (рентген, ультразвук, лазер и др.) и для любых типов материалов, если получен- 
ные данные представляют собой цифровое (растровое) изображение. 
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Ускоренный препроцессинг в задаче поиска подстрок в строке` 


А. В. Мазуренко', Н. В. Болдырихин”” 
'ООО «ДДОС-Гвард», г. Ростов-на-Дону, Российская Федерация 


* Донской государственный технический университет, г. Ростов-на-Дону, Российская Федерация 
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Введение. Бурное развитие таких систем, как Уапдех, 
Ооозе и пр., предопределило актуальность задачи поиска 
подстрок в строке. На сегодняшний день активно иссле- 
дуются подходы к ее решению. Эта задача используется 
при создании систем управления базами данных, поддер- 
живающих ассоциативный поиск. Кроме того, она приме- 
нима при решении вопросов информационной безопасно- 
сти, создании антивирусных программ. Алгоритмы поиска 
подстрок в строке используются в задачах обнаружения, 
основанного на сигнатурах. 

Материалы и методы. Решение задачи базируется на 
алгоритме Ахо — Корасик, который представляет собой 
классический способ осуществления поиска подстрок в 
строке. Вместе с тем применен новый подход в части, 
касающейся предварительной обработки. 

Результаты исследования. Показана возможность постро- 
ения функции перехода и суффиксных ссылок при помо- 
щи суффиксных массивов и специальных отображений. 
Исследована взаимосвязь между префиксным деревом и 
суффиксными массивами. Это дало возможность разрабо- 
тать принципиально новый способ построения функций 
перехода и ошибок. 

Полученные результаты позволяют существенно сокра- 
тить время, затрачиваемое на предвыборную обработку 
множества строк образцов при использовании целочис- 
ленного алфавита. 

В статье приведено восемь алгоритмов. Оценены разрабо- 
танные алгоритмы. Полученные результаты сопоставлены 
с ранее известными. Доказаны две теоремы и восемь 
лемм. Приведены два примера, иллюстрирующие особен- 
ности практического применения разработанной процеду- 
ры препроцессинга. 

Обсуждение и заключения. Предложенная в данной статье 
процедура препроцессинга основывается на связи между 
суффиксным массивом, созданным на основе множества 
строк образцов, и построением функций перехода и оши- 
бок на начальных этапах работы алгоритма Ахо — Корас- 
ик. Такой подход отличен от традиционного и требует 
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использования алгоритмов, позволяющих построить суф- 
фиксный массив за линейное время. Таким образом, опи- 
саны алгоритмы, позволяющие существенно сократить 
время на предварительную обработку множества строк 
образцов при условии использования определенного типа 
алфавита по сравнению с известным подходом, предло- 
женным А. Ахои М. Корасик. 

Результаты исследований, приведенные в статье, могут 
быть применены в антивирусных программах, использу- 
ющих поиск сигнатур вредоносных информационных 
объектов в памяти вычислительной системы. Кроме того, 
данный подход к решению задачи поиска подстроки в 
строке позволяет значительно ускорить работу систем 
управления баз данных, применяющих ассоциативный 
поиск. 


Ключевые слова: поиск подстроки, алгоритм Ахо — 
Корасик, префиксное дерево, суффиксный массив, поиск 
информации, функция ошибок, функция перехода. 
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Введение. В настоящее время существенно возрастает важность обеспечения кибербезопасности рас- 
пределенных информационных систем и отдельных вычислительных средств [1]. Диапазон таких задач доста- 
точно широк [1-10]. Особое место занимают вопросы создания эффективного антивирусного программного 
обеспечения (ПО). Одной из важнейших задач, решаемых таким ПО, является поиск подстрок в строке [1, 5, 6, 
10-13]. 

Материалы и методы. Задача поиска подстрок заключается в нахождении всех строк в тексте Т общей 
длиной т, совпадающих с каким-либо образцом из заданного множества образцов Р. Положим, сумма длин 
всех элементов Р, состоящих из символов алфавита /, равна п. Решение этой задачи было предложено А. Ахо и 


М. Корасик [6, 10]. В их алгоритме время предвыборной обработки составляет О(т! |, а время поиска — 


О(т! |+ К) . Здесь К — количество совпадений, найденных в тексте со строками, принадлежащими множеству 


образцов. 

В настоящее время задача поиска подстроки в строке активно исследуется по двум причинам: 

— поисковые системы активно развиваются [11]; 
— процесс обнаружения в антивирусных программных продуктах основывается на сигнатурах [1]. 

В этой связи были созданы алгоритмы, которые приходится выбирать в зависимости от определенных 
потребностей пользователя. Новейшие результаты, полученные при решении задачи поиска множества под- 
строк, описаны в работе [13]. 

Представленные в данной статье результаты основываются на связи между суффиксным массивом, со- 
зданным на основе множества строк образцов, и построением функций перехода и ошибок на начальных этапах 
работы алгоритма Ахо — Корасик. Такой подход отличен от традиционного и требует использования алгорит- 
мов, позволяющих построить суффиксный массив за линейное время. Итак, в статье описаны алгоритмы, с по- 


мощью которых время выполнения предвыборной обработки сокращается до О(п) 


* . 
Пусть задан алфавит [, множество образцов Р={В,Ь,...,В;| ‚ ге Ре/ , 1=ЬК. Обозначим через 
К < < 
п = > ИВ . Будем считать, что алфавит /[ представляет собой ограниченный промежуток целых чисел. Грани- 
1= 


* 
ца может зависеть от длины рассматриваемой строки 5Е/ либо представлять собой промежуток [0, с] ‚ гдес 


—щ натуральное число: с> |5 . Пусть & Е / — пустая строка. 


Пусть вою — функция перехода, а {аЙиге — функция ошибок. Приводимые модификации касаются ме- 
тодов построения упомянутых функций, используемых в алгоритме Ахо — Корасик [6, 10]. 
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* 
Положим, 5и/А77(5) — некоторый алгоритм построения суффиксного массива для строки 5Е/Г за 


линейное время. Описание таких алгоритмов можно найти, например, в [12—15]. 


* 

Пусть х,уЕГ . Тогда, [ср(х,у) — наибольший общий префикс строк хиу. 

Рассмотрим строку 5ЕГ у И $[5[0]5[1..5[и-—1]. Пусть $[5 [152 +1]...$[ Л| — подстрока 5, включаю- 
щая символы от? до], где 1< у, 1,1 =0,п-1. Обозначим через р, суффиксный массив, соответствующий 
строке 5. Пусть 

ру = Ру [ 50] ру... [п], 
то есть 5] 5[25[0]..5["-1] е 5[ 5 [2;[1]...5 [и -П| нь 5[ 5 [5 [п —1]-..5 [п -П| 

Для построения суффиксного массива будет использован алгоритм, описанный в [15]. 

Положим, а; #1, и *а,, 1<1<]<^+1, @ <9> <...< 91. Пусть УБеГ а; <Б, где 1<1<К+1. 
Допустим, что Р*@, арйа = {0,0 ,. бб} | 

Алгоритм обработки суффиксного массива р. 

Здесь зе Г: 5 =01 Аа»... Во, ВР; еГ, 1<1<. 

Ааарапоп (5, т арйа) 

|. пеу’_аттау < 8 


2.10 (: < арйа #++) { 





<} 





3.1 < 0 
4. ийЦе (5 [55 [11...59 -П [7 [2 арпа) { 


Е пе» аттау| [7] <—5| $р, [2]... 8] 1 
6. < 1+1 

7. } 

8. } 


9. огаегеа_ Пя [0] < пез? атгау [о] 


10. /ог (1—1; < [$] -|а/рйа #++) { 





11.7=—0 
2, И (пем’_ аттау[Й == пеу’ _ атау[!—!]) { 


13. от4еге4_1154 [1 | < пеу аттау | | 


14. < 7+1 

15. } 

16. } 

[7. гаити отаегеа_ [5 


Лемма 1. Пусть Р = |В,Р, | ‚ 5=09В0>Р...а, В, 0-1. Гогда алгоритм Адараноп строит массив 


лексикографически упорядоченных суффиксов образцов, принадлежащих Р, за время 0(|5 _ [рва] | 


Доказательство. В цикле 2—8 выполняется построение массива пеу’ а’тау, 1-й элемент которого пред- 
ставляет собой такой префикс соответствующего суффикса 5, который включает в себя все символы данного 
суффикса, начиная с нулевой позиции до его первого элемента, принадлежащего множеству арйа. При этом 


при помощи с иксного массива перебираются все с иксы 5 согласно их лексикографическому поряд- 
К) 


ку. Таким образом, массив пеу’ аггау состоит из всех суффиксов образцов, принадлежащих Р согласно их лек- 
сикографическому упорядочиванию, причем возможно повторение некоторых суффиксов. Заметим, что из рас- 
смотрения исключаются все строки, начинающиеся с символов, принадлежащих массиву а[рйа, то есть первые 


|@рйа суффиксы. Далее в цикле 10-16 при помощи массива пеу’ аггау строится массив огаегей_ 151 путем ис- 


ключения повторяющихся строк. Для этого в силу лексикографического порядка следования строк достаточно 
проверить, совпадает ли рассматриваемая строка с предыдущей. 
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Цикл 2-8 выполняется за время 0(|$]|-|а1рйа]) , поскольку из рассмотрения исключаются все строки, 
начинающиеся с символов, принадлежащих массиву арйа. В цикле 10—16 происходит |5] — |а/рйа] сравнений 


строк. Таким образом, получаем асимптотическую оценку времени работы алгоритма 0(|$]|-|айрйа]}. Лемма 


доказана. 
Алгоритм разбиения согласно лексикографическому упорядочению 
Здесь 5 — массив лексикографически упорядоченных строк. 


РапаС (5) 

1. зи [0] < 0 

2.1 <0 

3. Лог (1=—0;1<|$]-11++) { 


4. И (Я = [1+1] { 
5. зи6 [1 | < {+1 

6. < 7+1 

7. } 

8. 5иБ [7 — 5 

9. гаити 5иБ 


Лемма 2. Алгоритм ДапаС на основе массива лексикографически упорядоченных строк 5 строит массив 
5иб, состоящий из целых положительных чисел, указывающих индексы, соответствующие первым строкам сре- 


ди строк, у которых равны первые символы, за время 0(|) 


Доказательство. Граница, соответствующая первому символу, начинается с 0, что соответствует при- 
сваиванию, выполняемому на шаге 1. В цикле 3—7 последовательно сравниваются первые символы 1-Й и (1 + 1)- 





й строк, где 1=0, 5] —2. Если символы не равны, то в массив 5иБ записывается начало границы, соответствую- 


щей следующему символу. В противном случае выполнение цикла продолжается. Правая граница последнего 
символа соответствует количеству строк в массиве 5 (шаг 8). 


Сравнение на шаге 4 происходит за время (1) ‚ как и запись на шаге 5 и инкрементирование на ша- 


ге 6. Таким образом, цикл 3—7 выполняется за время 0(|5]) . Лемма доказана. 

Алгоритм построения первой связи 

Здесь "ее — дерево, [ех _ иогах Е [ ” ‚ Пик _ пит — номер некоторого символа строки [ех _ мог ‚ у — 
порядковый номер нового узла, который присоединяется к узлу с порядковым номером иоае _ питБег. 

ВийаР15Илик (тее&, [ех _ могау &, у&, Пик _ пит, поае _ питфег) 

1. пеу? тее.поае |] 

2. (тее.по4е [|] „5сие <— [ех имога$ [/ех _ иога$О]../ех _ иот[Итк] | 

3. пем’ тее.по4е [по4е _ питье | „Пик <— тее.поае |] 


4. (гее.по4е [по4е — питье | „Пик.5утбо[<— [ех_мота5 | ШК _ пит | 


5. У<У-| 
Лемма 3. Алгоритм ВийаЕ5Илик осуществляет построение нового узла с порядковым номером у и ду- 


ги, ведущей от узла пойе_питьег к новому узлу ъ, в дереве {тее за время (1) | 
Алгоритм построения связи с подстрокой 
* < < 
Здесь ее — дерево, [ех _ мот ЕГ ‚ у — порядковый номер нового узла, который присоединяется к 


узлу с порядковым номером 5а7. 


Вийа$иБ5 те Лик (тее&, [ех _ иог45&, у&х, 5ватт) 


1. юг (^ <— маи; К < [ех _ мот 





++) { 
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2. пе’ тее.поае |] 
3. тее.по4е |] „5 <— [ех_ уога$ [/ех _ иога$О]../ех _ мот] 
4. пеи’ гее.поае [у — || „Пик <— тее.по4е [] 


5. гее.по4е [у _ || „Пик.5утро[ < [ех_мота5 [К] 


6. у<у+1 
а. 
Лемма 4. Алгоритм Вийа$ирятте[алиК осуществляет построение новых узлов в дереве {тее, соответ- 


ствующих всем префиксам строки [ех_уота$, начиная с префикса [ех_ иога$ [/ех _ уого[О]..41ех _ ог [ат | 
за время О([ех _ мот — 51) 

Алгоритм построения последней связи 

Здесь "ее — дерево, [ех _ мога Е 1 " ‚ у — порядковый номер новой дуги, [ — алфавит. 

ВиИаГазИлик (тее&, [ех _ мога$, у, [ ) 

1. пе? (тее.поае [о] ИтЕ[у| — п'ее.поае [0] 

2. зутро[ < © 


3. юг (. — 0;1< [ех _ мот 





#++) { 
4. зутБо[5 [1] < ([ех_ мога$ [О] 
5.] < 0 

6. ог (= 0;1 < 





#++) { 

7. (ПЛ я зутБов) { 

8. п'ее.по4е [0] ИтЕ[у] 5утро [7] <= ПП 
9. 1<— 7+1 

10. } 

11.} 


Лемма 5. Алгоритм ВиИаГазИлик строит при корневом узле петлю. Ее пометка соответствует множе- 
ству символов, по которым невозможно перейти в другие узлы дерева ее из корневого узла за время 


О([ех мог] +| |) 


Алгоритм построения переходов 
Здесь [ех _ иог4$ — массив лексикографически упорядоченных строк. 


Стеше[лик (1 ех _ ига 


1. и < © 
2. виБ < ОапаС (1ех — мот) 


3. ух 1 
4. тее < © 


5. тее.по4е [о] „5(е < 5 

6. ог (#<— 0;1 < [ви] —11++) { 

7. ВиНаЕи$Илик (те, [ех _ иога5 [ий | ‚у, 0, 0) 

8. ВиНазибу те ик (те, [ех _ мот [и [| ‚У, 1) 

9. ог (1 < 5иИ ПБ] < зи +1; ++) { 

10. етр < Ср (1ех_ мое 7 — |,[ех _ мот +1 

11. < тее.зекиеМитьБег( [ср (1ех _ иога$[ 7 -—П[ех _ мога$[ 7] ) 
12. ВиПаЕи$Илик (1тее, [ех _ мога$[ 7], у етр, 2) 
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13. Вий$ир5 те Лик (ее, [ех _ иога5 | 7 | У, етр} 

14. } 

р. 

16. ВипаГа;Илик (те, [ех _ мот, у[ех _ иота 

[7. гаиги 1тее 

Лемма 6. Алгоритм СтежшеМик осуществляет построение префиксного дерева {ее с петлей при корне- 


|[ех в ига -1 


ый |[ех _ нотам 


вом узле за время с 


Доказательство. На шаге 2 выполняется алгоритм ДапаС (см. лемму 2), после чего на шаге 5 создается 
корневой узел с порядковым номером 0 дерева {тее. Его состояние полагается равным пустой строке &. Рас- 
смотрим цикл 6-15 на 1-м шаге работы. 

На шаге 7 при помощи алгоритма ВийаР 5 Илик создается узел, состояние которого соответствуют пер- 


вому символу строки [ех _ ога [и [7 |. Учитывая построение массива 5иб, можно утверждать, что такой сим- 


вол прежде не встречался среди первых символов предыдущих строк. Тогда на шаге 8 реализация алгоритма 
ВиПа5ибуттеЛик последовательно создает узлы, состояние которых соответствует всем префиксам строки 


[ех _ мог [и] ‚ исключая префикс, построенный в предшествующей строке. 
В цикле 9-14 при помощи алгоритмов ВийаЕ $ Илик и ВиНа5ибятте[атК выполняем те же действия со 
строками, лежащими в целочисленном промежутке [би -1, 5иБЁ +1-1]. Поскольку у каждой такой строки 


есть общий ненулевой префикс с предыдущей строкой, то алгоритм немедленно переходит в состояние, соот- 
ветствующее наибольшему общему префиксу, начиная с которого необходимо строить новые узлы. На шаге 16 
при помощи алгоритма ВийаГазИлик создается петля при корневом узле. 

Шаги 12 и 13 выполняются за время 


О(1)+ О(ех_ у’отеб И _ р (1ех_ у’отеб [7 _ 1] „[ех _ мог [1 _ у = 


= О(1ех_ мот] — Иер (1ех _ уогаб [7 — 1] „[ех _ иога$ [1] 


Таким образом, из лемм 2, 3 и 4 следует, что цикл 9-14 выполняется за время 


Ы 
ое иоаз[ 7 = «р ( (1ех _ мог [= 1], [ех _ нота Л 


Цикл 6-14 выполняется за время 


с[У. С [ех _ иогаб БынА] + 


ое иотаз [| = «р [ех _ мой [7 — 1], [ех _ нон} 
= 1=5иб 


[би] —1 
=оО [ех_ мото [ зи [7] |--|ех_ мое] |. 


1=0 


|ех м ига 
Из леммы 5 следует, что шаг 16 выполняется за время ох вонь ) 
1= 


[ех _ мога$ 0] 


Таким образом, получаем асимптотическую оценку времени работы алгоритма 


О([ех — этот} + с (У. ка [ех _ иогаб [5и ] — [ех — нота Е 
ра 


+О [е — иоте] + >. т в [ех _ ОГ 0] =О [> в ры [ех _ Ога 0] 
1- [= 


Лемма доказана. 
Алгоритм построения функции 200 
Здесь Р — множество строк образцов. 


СопятгисЮою (Р) 


1. 5 < 0 Ао>Р>...04. В; 041 
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2. р; < зиЙАи(5) 

3. арйа <— {0,92 .. об! 

4. отаеге4_П51 < Ааарайоп (5, вы рва) 
5. < 0 

6. [ех_ ога < © 


7. Р [епей < В ь 








} 


8. Лог (#=0;1 < 





++) 
9. ((огегеа _ И еЕР_[епо1й) апа (отаегеа _ ПУ] Е Р)) 
10. [ех_иога$ [71++] < огаегей _П5П 


11. гою < Сгеше“лиК(ех _ иога$) 


12. гешти 2010 


[4 
Напомним, что и = > [В для множества образцов Р= {В, Еее В, } 
1= 


Теорема 1. Алгоритм Сопягисою строит функцию го за время О(п) 


Доказательство. На шаге 2 строится суффиксный массив р; для строки 5. На шаге 4 при помощи алго- 


ритма Аааранйоп в массив огаегей_ [5 записываются все суффиксы строк, принадлежащих множеству образцов 
Р. При этом не исключены повторения. В цикле 8—10 строится массив [ех_мога$, содержащий суффиксы, при- 
надлежащие Р, расположенные в лексикографическом порядке без повторений. На шаге 11 происходит постро- 
ение префиксного дерева с петлей при корневом узле на основе строк, которые содержатся в массиве [ех_и’огаб. 
Структура данных, возвращаемая алгоритмом Стеше[лт^, в точности задает функцию 20. 


Шаг 2 выполняется за время О(п+К+1) [12]. Из леммы 1 следует, что шаг 4 выполняется за время 


О(п ++1-А-— 1) = О(п) . В цикле 8—10 рассматриваются только строки, длина которых равна длине какого- 


либо образца. 


Таким образом, потребуется не более О(п) проверок для нахождения образцов Р. Из леммы 6 следует, 


|[ех № мот -1 


а |[ех _ нота =О(п). Так как А < п ‚ получаем асимптоти- 
1= 


что шаг 11 выполняется за время с] > 


ческую оценку времени работы алгоритма О(п) + О(п +А+ 1) = О(п) . Теорема доказана. 
Результаты исследования 
Пример 1. 
Пусть Р= { опе, оп, опсе, сей, иИ,еуе, пеаг} .‚ Тогда 
5 = ал опеаропозопсеодсеПаиПовеуеал пеагов . (1) 
В табл. 1 представлен результат работы алгоритма построения функции го при поступлении на вход 


строки 5 (1). 
Таблица 1 


Структура префиксного дерева 


Номер Состояния, в которые ведут ветви Символы на ветвях Дик от 
Состояние ипоае 
поае Дик от поае поае 


1.с;:2.е;: 3. 1:4.п;5.0 1.с;:2.е;: 3. 1:4.п;5.0 
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Номер Состояния, в которые ведут ветви Символы на ветвях Дик от 
Состояние ипоае 
поае Дик от поае поае 


и 1. 1.1 
0 [7 1. м 1.1 
т 
2 1.е 
3 1. пеа 1.а 
4 1. пеаг 1.7 
5 
6 1. оп 1. и 


й 1. опс; 2. опе 1.с;2.е 


1. оисе 1.е 


Пусть 9 =а; В. ах... Во В) — зеркальное отражение строки 5. 





Алгоритм построения функции ошибок /аЙиге 
Здесь Р — множество строк образцов. 


Бабеби(Р) 

1. < а; 1В. ау... Ра В) 

2. р; < эиЙ Аг ($) 

3. арйа < {0,9 . | Зорб 

4. огдеге4_ 11$: <— Адарайоп (5 ‚р ‚рйа) 
5. Дик < © 

6. ог (<= 0;1<|5 
7. апйк[й < 8 

8. зи < ДапаС (ог4еге4 — [51 
9. 5" <— © 

10. о" (1—0; < [би —11++) 





++) 


И. Юг (1 < зи Д; 1 < зи + - В ++) 

12. 5" [7] <= /ср(огаеге4 _ [58 7], огаеге4 _ [5 + | 
13. ог (; < 0;1< [5и | -11++) 

14. ог (К < 5и +|-БК > зи ПК --) { 

15. ог (1<5и МИ; <; ++) 

16. тт_еетет < т (5-1, 5—2]... Л) 
И.И (пит _ е/етеи! = |огеге4 _ 15 7) 

18. шш етр [7 —5и Й] <— тт_еетеп! 


19. } 


20. тах_еетепи <— тлах (ти _ {етр[0], пи _(етр[],...,пиа етр[м]) 
21. найти тах_таех: шш_етр [лах_иа4ех | = тах @етет 
22. [тк [К] < огаегеа _ [5 [лах_ таех + 5и | 


23. } 
24. дог (<= 0;1 < [тик 
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25. Пи ИО] < огаегеа _ Пя |` ; зеркальное отражение строки 


26. ПИ] < тым ; //зеркальное отражение строки 
Ту 


28. гаитп Пик 

Замечание. В строке 20 < 5и 1+ П-5и -1. 

Теорема 2. Алгоритм Рабези} строит функцию ошибок /аЙиге за время О(и). Доказательство. На ша- 
ге 1 строим массив символов, содержащий зеркальные отображения строк, принадлежащих множеству образ- 
цов Р, и некоторые уникальные символы. На шаге 2 строим суффиксный массив р; для строки 5. На шаге 4 
при помощи алгоритма Адарайпоп в массив огаегеа_ [5 записываются все суффиксы строк, принадлежащих 
множеству образцов Р (множество образцов, состоящее из зеркально отраженных строк Р), при этом не ис- 
ключены повторения. На шаге 8 выполняется алгоритм РапаС (см. лемму 2), после чего в цикле 10—12 находим 
длину наибольшего общего префикса между строками, у которых совпадает первый символ. Записываем ре- 
зультат в массив 51". Заметим, что эта величина равна нулю у строк, для которых данное условие не выполняет- 
ся. В цикле 13—23 строится специальное отображение между строками, у которых совпадает первый символ. 
Опишем это отображение. Зафиксируем некоторую строку 


5 е огаегей _ [51 | ог4еге4 _ [5 [ий | .... огаегеа _ 151 [5 П + - || 
Рассмотрим множество, состоящее ИЗ строк, принадлежащих 


огдеге4 _ 151 | опдеге4 _ 151 [ий |,..., ог4еге4 _ 15 [56 1+1-1]. Их длина равна длине наибольшего общего 


префикса с 5, исключая само 5. Из этого множества найдем строку 5’, которая имеет наибольшую длину, и по- 
ставим ее в соответствие 5. Очевидно, что построенное отображение является биекцией при условии 5’=8. 
Полученный результат записывается в массив т/лик. В цикле 24—27 при помощи массива т/лик явно указываем 
построенное отображение, при этом зеркально отражая каждую из строк. Таким образом, ставим узлу $ пре- 
фиксного дерева, построенного на основе массива образцов Р, в соответствие узел 5’. Его состояние равно 
наибольшему собственному суффиксу 5 , который встречается среди множества состояний рассматриваемого 
префиксного дерева. Но согласно определению функции ошибок /аЙиге это и есть искомый результат. 


К 
Пусть и = > [В . Шаг 2 выполняется за время О(п+к-1) [12]. Из леммы 1 следует, что шаг 4 вы- 
1 = 


полняется за время О(п +А+1-—А^ —1) = О(п) . Цикл 6-7 выполняется за время О([5 | = О(п) . Из леммы 2 следу- 


ет, что  шаг8 выполняется за время О(п) . Цикл 10-12 выполняется за время 
Ь-2 Ы 1-1 
[у Хиты — О(к ), М7 = 1. Цикл 14—23 выполняется за время 


5иб [1+1 к 5иб [1+ -5и 61-2 
О ‚ |= О ‚ |= О(5иБТ+П-5иЙ-Т), М] И 
их, ип" У 7—0 ',) (5 р - ) 1!) 


Тогда цикл 13—23 выполняется за время 
[5и6|--2 , , 
с[ У." (ви +0-змий-1) = О(5иь[[5и6| - 1) = (К). 
Поскольку тт <п, то цикл 24—27 выполняется за время О([тЬт]) = О(п). Итак, учитывая, что 


К <пи, получаем асимптотическую оценку времени работы алгоритма О(п) + О(п +А+ 1) + О(К) = О(п) . Тео- 


рема доказана. 
Пример 2. 
Пусть, как и в примере 1, Р= { опе, оп, опсе, сей, иИ, еуе, пеаг} . Тогда 
5 = озтаеполеуеав Иша Песадесповзпобрепос\ . (2) 
В табл. 2 представлен результат работы алгоритма построения функции /аЙиге при поступлении на 
вход строки $ (2). 
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Таблица 2 


Ложные связи между узлами 


ве в бе | ве | № ое. 


Для всех узлов, для которых на рис. 1 не указаны ложные связи, полагаем, что ложная связь ведет к 





корневому узлу. 
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Рис. 1. Префиксное дерево с ложными связями 


Обсуждение и заключения. Описана новая процедура препроцессинга в алгоритме Ахо — Корасик. 


Она выполняется за линейное время О(п). Исследована связь между суффиксными массивами и префиксным 


деревом, что позволило предложить иной способ построения функций перехода и ошибок. Полученные резуль- 
таты позволяют сократить время на предвыборную обработку множества строк образцов при использовании 
целочисленного алфавита. 
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